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1.  INTRODUCTION 


n*  • 


I 

Continuum  mechanics  codes  have  been  applied  to  a wide 
variety  of  impact  situations,  including  those  in  which  sevete 
material  damage  occurs  to  both  the  penetrator  and' the  target 
material.  It  is  important  to  include  in  such  calculations  an 
appropriate  description  of  the  material  damage  generation 
process. 

Recently,  rather  detailed  material  failure  models 
have  been  developed  at  SRI  International which  relate 
material  failure  to  the  nucleation  and  growtn  of  microscopic 
holes  or  cracks.  These  models,  known  collectively  as  NAG-FRAG 
(Nucleation  ^d  Growth-FRAGroentation) , were  developed  for  both 
ductile  materials  (the  DfkACT  model)  and  brittle  materials 
(the  BFRACT  model) . These  models  have  been  incorporated  into 
the  1975  documented  version  of  HELP'^'  to  produce  the  present 
NAG-FRAG  version  of  the  code  which  gives  a much  more  detr  ^ led 
description  of  material  failure. 

The  major  modification  to  the  sHELP  code  for  the  incorpo- 
ration of  these  models  involves  the  establishment  and  subsequent 
manipulation  of  "Lagrangian  damage  regions."  These  damage 
regions  encompass  pockets  of  failed  material  and  move  across 
the  Eulerian  mesh  with  the  local  velocity  field.  A damage 
region  is  formed  when  material  in  an  Eulerian  cell  first 
satisfies  the  failure  criterion. 

The  two  NAG-FRAG  subroutines,  BFRACT  and  DFRACT,  communi- 
cate with  the  external  HELP  code  by  calculating  the  stress  state 
of  the  material  in  the  damage  regions.  In  turn,  the  external 
code  provides  BFRACT  and  DFRACT  with  such  quantities  as  material 
density,  internal  energy  and  the  strain  increments  to  be  used 
in  the  calculation  of  the  damage  region  stresses.  Furthermore, 
unlike  the  properties  of  neighboring  Eulerian  cells,  the 
characteristics  of  contiguous  damage  regions  are  not  diffused 
because  there  is  no  material  convected  from  one  to  the  other. 

This  repor.  documents  the  changes  that  were  necessary 
to  Incorporate  the  SRI  NAG-FRAG  models  into  the  HELP  code. 

In  Section  2,  these  models  are  briefly  described,  and  in 
Section  3 the  modifications  that  were  made  to  the  HELP  code 
for  their  incorporation  are  discussed.  Section  4 discusses 
several  auxiliary  improvements  that  were  made  to  the  HELP 
code  concurrently  with  the  incorporation  of  the  NAG-FRAG  models. 
In  Sectiori  5 results  from  several  calculations  performed  using 
the  version  of  HELP  that  contains  the  NAG-FRAG  models  are 
presented.  These  results  are  compared  with  Lagrangian  calcu- 
lations performed  by  SRI.  Appendices  A and  B are,  respectively, 
a user's  guide  to  the  revised  code  which  supplements  the  1975 
documentation,  and  a new  set  of  input  forms. 
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By  incorporating  tha  NAG-PRAG  failtira  models  into  the 
Euler ian  HELP  code,  an  inf}roved  n\»Mirical  tool  has  been 
developed  for  the  simulation  of  the  armor  penetration  process 
which  involves  both  large  material  deformations  and  extensive 
material  failure.  This  numerical  tool  can  be  used  by  designers 
of  both  penetrator  and  anaor  systems  and  can  provide  detailed 
data  to  vulnerability  and  survivability  analysts. 


t 


2.  A DESCRIPTION  OF  THE  NAG-FRAG  MODELS 


The  NAG-FRAG  (Nucleation  Growth-FRAGmentation) 
models  incorporated  into  the  HELP  cocTe  were  developed  by  SRI 
for  predicting  the  damage  incurred  by  ductile  and  brittle 
materials  when  subjected  to  tensile  stresses.  The  damage  is 
assumed  to  result  from  the  generation  of  small  spherical  voids 
(ductile  failure)  or  disc  shaped  cracks  (brittle  failure)  and 
is  reflected  by  a degradation  of  material  strength  properties. 

Reference  [1]  describes  these  models  (DFRACT  for  ductile 
materials  and  BFRACT  for  brittle  materials)  in  detail.  The 
following  is  a brief  description  of  those  aspects  of  the 
models  for  which  an  understanding  was  necessary  in  the  develop- 
ment of  the  NAG-FRAG  version  of  HELP.  In  the  course  of  this 
description*  reference  will  be  made  to  the  TSR  array  which  con- 
tains those  material  properties  unique  to  the  NAG-FRAG  models. 


2.1 


VOID  NUCLEATION  AND  GROWTH  IN  DUCTILE  J4ATERIALS 


The  spherical  voids  responsible  for  damage  in  ductile 
materials  are  assumed  to  (1)  be  generated  at  a prescribed 
rate  if  the  tensile  pressure  exceeds  a nucleation  threshold 
and  (2)  grow  in  size  at  a prescribed  rate  if  the  tensile 
pressure  is  above  a growth  threshold.  During  nucleation* 
the  voids  are  assumed  to  be  generated  according  to  a given 
size  distribution. 

More  specifically*  the  nucleation  rate  N (number  of  voids 
generated  per  unit  volume  per  unit  time)  is  given  by  the  follow- 
ing relations  wherein  pressure  is  defined  to  be  negative  in  tension 

P-P 

no 


N = 


e 

o 


if  P < P 


no 


if  P > P 


no 


where 


P ■ tensile  pressure 

Pro  “ TSR(m,5)  = nucleation  threshold  pressure 
Pj^  = TS.A(m,6)  =»  nucleation  sensitivity  factor 

= TSR(m*4)  = nucleation  rate  at  threshold  pressure. 

The  first  index  of  the  TSR  array,  ra*  identifies  the  particular 
material  being  modeled. 
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The  void  radii  at  the  time  of  nucleation  are  assumed 
distributed  according  to  ' 


f (R) 


1 - 


-R/R„ 


where  f(R)  is  the  fraction  of  voids  having  a radius  less  than 
R and 


R^  * TSR(mi3)  > nucleation  void  radius  distribution 

parameter. 

Thus,  the  volume  of  voids  nucleated  per  unit  voltune  of  material 
during  a time  step  dt  is  (assuming  P remains  constant  during 
the  time  step) 


AV 


n 


At 


df  (R) 

"Ir 


dR 


Btt  N R**  At 
n 


The  void  radii  are  assumed  to  grow  according  to  the  law 
fP-P. 


4ti 


if  P < P 


go 


if  P > P ^ 
- go 


where 


go 


material  viscosity 

TSR(m,2)  * growth  threshold  pressure 


Thus,  if  a particular  void  has  a radius  R at  the  beginning 
of  the  time  step,  its  radius  at  the  end  or  the  time  step 
will  be  p_p 

''  -7^  At  if  P < P„„ 

4ti  go 


«o  « 


if  P > P„^ 
- go 
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As  a consequence  of  this  growth  law»  the  void  content  of  a 
mass  of  damaged  material  remains  unchanged  during  recompression. 
Given  a void  volume,  V = 4/3  n R^,  at  the  beginning  of  the 
time  step,  the  void  volume  at  the  end  of  the  time  step  will 
be 


V » e 
o 


Tn  “'■'’go 


The  factor  3/4n  is  input 


)At 

to  the  model  as  TSR(m,l). 


2.2  CRACK  NUCLEATICN  AND  GROWTH  IN  BRITTLE  MATERIALS 

The  nucleation  and  growth  of  cracks  in  brittle  materials 
is  treated  in  a fashion  much  like,  but  more  complex  than,  the 
voids  in  ductile  materials.  The  number  of  cracks  nucleated 
with  a particular  angular  orientation  is  assumed  governed 
by  a nucleation  rate 

°no 

°1 

e 


o 


where 

o,.  * stress  normal  to  the  plane  of  the  crack 

4"!' 

o * TSR(m,5)  ■ nucleation  threshold  stress 
1^0 

5 TSR{m,6)  = nucleation  sensitivity  factor 
N = TSR(m,4)  ■ nucleation  rate  at  threshold  stress. 


“ ‘ °no 


if  o, . > 

no 


In  this  section  of  the  report  and  in  the  NAG-FRAG  subroutines, 
both  the  pressure  and  the  total  and  deviatoric  stresses  are 
defined  to  be  positive  in  compression.  So  as  not  to  consider 
every  possible  angular  orientation,  five  bin  angles  are  con- 
structed as  shown  in  Figure  2.1.  For  exeunple,  all  cracks 
whose  faces  are  iiormal  to  the  direction  given  by 


45®  i 4)  i 135®  and  -22.5®  i (|>  i 22.5® 
or 

0®  < \|)  < 45®  and  0®  < ((>  < 360® 
are  included  in  bins  1 and  5,  respectively. 


Crack  radii  growth  ia  assumed  governed  by  the  equation 


1' 


' ' if  0 ” ’<  a L *- visas  s *■  -i  » l- 

♦♦r,  . 9®  6 liirr-nootf-j 


if  O.  . i. 


, I 


or,  analogous  to  void  growth. 


*^o  ® 


At 


if  o- < o„- 

♦'I'  go 


“ - V 


where 


jiri  1 

1 r-  1 
Y il.*ra>  ■ *B 


a « TSR(m,2)  « growth  threshold  stress, 
go 

The  cracks  are  assumed  to  open  elastically  to  form  ellipsoids 
having  a thickness  26  where  6 is  given  by 

-,1>C  , ’ • ' ‘i 

. . 4(1-v2)  * ' ■ - 

5 s R 


and 


' , t 


E > Young's  modulus 

> 

V K Poisson's  ratio 
- "»>'<-%'<»  • 

The  volume  of  a crack  is  therefore 


w ^ 1.2  . 16  ,l-v'.„3 

V - ^ R ^ - -y,  (-g — )R  0^^ 


J ** 


J.  D 


and  the  crack  volume  associated  with  bin  i is 


/ 
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where  as  before 


f(R)  ■ 1-e 


-R/R. 


Thus  It  is  assumed  that  a crack  radius  R remains  constant 
during  recompression  but  that  the  thickness  25  rnd  therefore 
the  volume  become  zero. 


2.3  DEGRADATION  OF  MATERIAL  PROPERTIES,  FRAGMENTATION 

The  nucleatlon  and  growth  laws  provide  an  estimate  of 
a void  volume  fraction  Vy.  The  shear  modulus  G and  the  yield 
strength  Y are  assumed  related  to  the  values  Gq  and  Yq  for 
undamaged  material  by  the  relation 


G - G^ 

Y » Y^ 
o 


where  the  factor 

* 15  (^) 

is  defined  by  a data  statement  in  the  models  to  have  the  value 
1.88  corresponding  to  a Poisson's  ratio,  v,  of  1/3. 

In  the  brittle  failure  model,  the' density  of  cracks 
together  with  a consideration  of  their  orientation  le/.ds  to 
an  estimate  of  the  extent  to  which  fragmentation  occurs.  This 
estimate  Includes  both  the  number  of  fragments  formed  and  their 
size.  As  fragmentation  becomes  conqplete  (the  quantity  FU2D 
goes  to  zero) , th;  .naterlal  can  no  longer  sustain  tensile 
pressures.  Fragmentation  is  governed  by  the  input  values  read 
into  TSR(m,10)  through  TSR(m,13). 


1 - p F 
1 - 4 V p 


2.4  EQUATION  OF  STATE 

It  is  assumed  that  the  equation  of  state  can  be  ade- 
quately represented  during  fracture  (because  of  the  small 
solid  strains)  by 

Pg  « Ap  + r pg  E 
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where 


\J‘ 

P > pressure  of  the  solid  materiel  (i.e.,  exclusive 
of  voids) 


A •>  material  bulk  modulus 
r * Gruneieen's  ratio 


'■-'tul  ■ 


^ y ■ •^s^^o  " ^s  density  of  the  solid 

material  and”p^  the  initial  density. 


Both  of  the  NAG-FRAG  models  determine  simultaneously,  during 
a time  step,  the  damage  quantities  (number  of  voids  or  cracks 
and  the  void  volume  fraction)  and  the  stresses.  This  is  done 
in  a con^lex  iteration  procedure  in  which  the  new  pressure 
estimate  must  be  carefully  initialized.  A detailed  discussion 
of  this  procedure  is  given  in  Reference  [1] . 
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3.  INCORPORATION  OF  THE  NAG-FRAG  FAILURE  MODELS 
INTO  THE  EULERIAN  HELP  CODE 


As  described  in  the  foregoing  section,  the  NAG-FRAG 
failure  models  provide  considerable  detail  as  to  the  extent 
and  type  of  damage  experienced  by  a given  volume  of  material 
over  a given  period  of  time.  In  a Lagrangian  code,  where  the 
computational  cells  are  associated  with  the  same  mas?  through- 
out the  calculation,  the  time-dependent  damage  characteristics 
(e.g.,  void  (or  crack)  number  and  volume)  of  each  cell  are 
easily  computed.  On  the  other  hand,  in  an  Eulerian  code, 
such  as  HELP,  where  the  grid  remains  fixed  and  the  mass  is 
transported  from  one  computational  cell  to  another,  these 
time-dependent  damage  characteristics  become  severely  dif- 
fused if  they  are  transported  as  the  momentum  and  energy 
quantities  are  transported. 

To  avoid  this  diffusion,  damage  regions,  which  deform 
with  the  mass  in  the  manner  of  a Lagrangian  cell,  have  been 
incorporated  into  the  HELP  code.  The  stresses  and  damage 
characteristics  of  these  d2unage  regions  are  computed  by  the 
NAG-FRAG  subroutines,  DFRACT  and  BFRACT.  The  Eulerian  cells 
that  are  completely  contained  within  the  boundaries  of  a 
single  damage  region  assume  the  strers  state  of  that  region; 
whereas,  those  cells  con^letely  outside  the  damage  region 
retain  the  stress  state  assigned  to  them  by  the  equation  of 
state  for  the  undamaged  material.  A cell  intersected  by  a 
damage  region  boundary  is  assigned  a stress  which  is  a volume 
weighted  average  of  the  stresses  associated  with  damaged  and 
undamaged  portions  of  the  cell. 

Once  the  stresses  of  the  Eulerian  cells  are  so  deter- 
mined, the  equations  of  motion  and  energy  and  the  transport 
of  mass,  momentum  and  energy  are  executed  as  before. 

To  insui'e  the  success  of  this  approach,  several  refine- 
ments were  made  to  the  computations  of  transport  volumes  and 
of  deviatoric  strain  and  stress  increments  in  HELP.  These 
refinements  along  with  a description  of  the  algorithms  for 
generating  ani  vipdating  the  damage  regions  are  discussed  in 
the  sections  th  it  follow.  Additional  information  of  the  type 
required  by  a u>er  of  the  code  is  given  in  Appendix  A. 


3.1  CREATION  AND  PROPAGATION  OF  DAMAGE  REGIONS 

A damage  region  is  created  when  material  in  an  Eulerian 
cell  first  fails,  i.e.,  the  cell's  pressure  (DFRACT)  or 
stresses  normal  to  bin  angles  (BFRACT)  are  sufficiently 


tensile  to  nucleate  voids  or  cracks.  Special  damage  region 
tracer  particles  are  placed  around  the  newly  failed  material 
and  thereafter  move  with  the  local  velocity  field.  These 
particles,  in  so  far  as  they  represent  the  boundary  of  this 
failed  material,  provide  a means  of  tracking  the  material 
as  well  as  measuring  its  volumetric  and  shearing  strains  as 
it  continues  to  deform. 

{ 


3.1.1  Movement  of  Deunage  Region  Tracer  Particles 


The  local  velocity  field  of  a damage  region  tracer 
particle  is  defined  in  the  same  manner  as  that  for  a material 
interface  tracer  particle.  The  velocity  components  of  those 
Eulerian  cells  overlayed  by  a rectangle  of  cell  dimensions 
that  is  centered  about  a given  tracer  particle  are  averaged 
to  give  the  velocity  components  of  the  given  tracer  particle. 
At  most,  four  Eulerian  cells  are  overlayed.  The  cosponents 
of  the  i^^  cell  are  weighted  by  the  area  common  to  it  and 
the  rectangle 


u 


lUiA. 

E A, 


E V. A, 


ta: 


This  procedure  gives  a spatially  continuous  velocity  field 
for  tracer  particle  motion. 


3.1.2  Creation  of  Damage  Regions  for  Pure  Eulerian  Cells 

If  the  *'ulerian  cell  which  has  just  failed  is  pure, 
i.e.,  contains  only  one,  undamaged  material,  the  boundary  of 
the  new  damage  region  coincides  with  the  boundary  of  the 
entire  cell.  Five  damage  region  tracer  particles  are  placed 
at  the  corners  of  the  cell  in  a counter-clockwise  order  with 
the  fifth  particle  being  a repeat  of  the  first. 
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3.1.3  Creation  of  Damage  Regiona  from  Multlicaterial 
^'lerlan  Celia 

When  one  of  the  materials  in  a multiinaterial  Eulerian 
cell  first  satisfies  the  failure  condition,  the  boundary  of 
the  new  damage  region  coincideu  with  the  boundary  of  that 
material  in  the  cell.  Damage  region  tracer  particles  are 
placed  at  the  intercepts  of  the  material  interface  with  the 
cell  boundary,  at  the  Included  corners  of  the  cell,  if  any, 
and  along  t^e  material  interface,  coinciding  with  the  material 
interface  tracer  particles.  The  last  particle  closes  off 
the  region  by  being  a repeat  of  the  first. 

t ^ 


• Damage  region  tracer 
particles  ' w-  ■ 

X Material  interface 
tracer  particles 

Material 

Interface 


3.1.4  Enlargement  and  Reduction  of  Existing  Damage  Regions 

A partially  damaged  cell  is  one  which  contains  the 
boundary  of  one  or  more  damage  regions  and  therefore  con- 
tains both  damaged  and  undamaged  materials.  During  each  time 
step  the  stress  state  of  the  undamaged  material  is  independently 
assessed  and  the  appropriate  failure  criteria  is  applied  to 
it.  When  the  failure  criterion  is  satisfied,  one  of  tx^o  courses 
is  taken:  either  the  newly  damaged  material  is  included  in 
an  existing  damage  region,  or  a new  d2unage  region  is  created  ^ 
to  include  the  previously  as  well  as  the  newly  failed  material; 
the  damage  regions  associated  with  the  previously  damaged 
material  are  reduced  accordingly.  The  first  option  is  chosen 
if  the  newly  damaged  material  occupies  less  than  half  the  cell 
volume,  or  in  the  case  of  a multimaterial  cell,  less  than  half 
the  volume  associated  with  the  given  material;  otherwise,  the 
second  option  is  chosen.  Figure  3.1  illustrates  these  two 
occurrences  for  the  case  when  the  partially  damaged  cell  con- 
tains only  one  material.  Note  that  in  the  case  where  an  existing 
damage  region  intersects  the  cell,  the  one  which  has  the  largest 
part  of  the  cell  within  its  boundaries  is  the  one  which  is 
enlarged;  the  others  are  reduced. 


Case  1:  A pure  Eulsrlan  cell  with  lees  then  one-helf  of  its 
material  previously  damagM  • a new  damage  region 
is  created,  and  the  overlapping  damaged  regions  are 
reduced . 


Case  2:  A pure  Eulerlan  cell  with  more  than  one-half  of  its 
material  previously  damaged  - one  of  the  existing 
damage  regions  is  enlarged,  and  the  ot.'terCs)  are 
reduced , 


Figure  3.1.  The  creation  or  enlargement  of  Lagrangian  damage 
regions. 


To  further  illustrate  the  way  in  whit  1 damage  regions 
are  created,  enlarged  and  reduced,  Figures  : .2  and  3.3  repre- 
sent a portion  of  an  Eulerian  grid  i r which  all  ol  these  op>l:ioiis 
are  exercised  and  the  damage  occtirs  in  inult:.!naterial  cells. 
Initially,  three  deunage  regions  exist,  as  sNown  in  Figure  32. 
The  next  figure  shows  the  configurar.icn  of  damage  regions 
which  is  produced  ^en  both  materials  in  all  of  the  mult. viu iter j.tl 
cells  fail  simultaneously  (an  unlikely  occurrence,  presented 
here  only  to  illustrate  the  various  mechanisms  for  generating 
damage  regions)  The  following  lists  the  steps  taken  to  pro- 
duce the  configuration  in  Figure  3.3: 

1.  Region  #1  was  reduced  to  exclude  cells  (2,1),  (2,2) 
and  (1,2)  since  less  than  half  of  Material  A in  any 
of  those  cells  was  previously  damaged. 

2.  Region  #2  was  enlarged  to  include  all  of  material 
B in  cell  (2,2)  since  it  had  already  encompassed 
more  than  half  of  material. B in  that  cell. 

i '1 

3.  Region  #3  was  eliminated  altogether  when  it  was 
reduced  twice;  once  when  region  #4  was  created 
and  again  when  region  #6  v t created. 

4.  Regions  44  through  #12  are  now  regions;  however, 

#4,  8,  10  and  6 encompass  some  material  that  was 
damagod  previously  and  belonged  to  regions  fl  or 
#3.  Regions  15,  7,  9,  11  and  12  encompass  only 
newly  failed  material. 

This  example  also  illustrates  the  motivation  for  en- 
larging and  reducing  existing  damage  regions  instead  of  always 
creating  new  ones.  Such  a procedure  minimizes  the  number  of 
damage  regions  that  are  much  smaller  than  a cell;  although  as 
regions  #5  and  #9  exemplify,  the  creation  of  small  regions 
is  not  always  avoided.  In  addition,  this  procedure  minimizes 
the  formation  of  damage  regions  having  very  irregular  shapes, 
such  as  the  undamaged  portion  of  material  A in  cells  (1,2)  ai.d 
(2,2)  in  Figure  3.2. 

The  boundary  of  an  existing  damage  region  is  enlarged 
or  reduced  by  removing  tracer  particles  interior  to  the  newly 
failed  material  as  well  as  by  adding  tracer  particles  at 
appropriate  cell  corners  and  at  intercepts  of  the  damage  region 
and  the  cell  boundaries. 

When  a region  is  reduced,  its  mass  >s  adjusted  as 
follows: 

« M?  - p.  AV 
1 1 1 
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Figure  3.3.  Configuration  of  damage  regions  after  both 
materials  fail  simultaneously  in  all  of  the 
multiroaterial  cells.  ^ 

1 - ' , . 
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where  p.  is  the  current  maee  density  of  the  region  and  AV 
is  the  ^lume  by  which  it  is  being  reduced.  Likewise,  when 
a region  is  enlarged,  its  mass  becomes  ^ 

- mJ  + pj  AV 

where  p . is  the  current  mass  density  of  the  newly  failed 
material  and  AV  is  its  volume. 

Enlarging  a dtunage  region  requires  all  of  the  material 
properties,  not  just  the  mass,  to  be  redefined.  The  damage 
parameters  and  the  stress  conqponents  of  the  enlarged  region 
are  redefined  to  be  a volume  weighted  average  of  those  quanti 
ties  associated  vrlth  the  region  before  it  was  enlarged  and 
those  quantities  associated  with  the  material  that  has  just 
failed  and  just  been  added  to  the  region.  For  example,  the 
pressure  of  the  enlarged  region  becomes 


- (P"  + P”  AV)/(V^  ♦ AV) 

where  1 and  j subscripts  refer  to  the  existing  and  newly 
damaged  materials,  respectively.  On  the  other  hand,  the 
specific  internal  energy  of  the  enlarged  region  is  a mass 
weighted  average  of  the  previously  and  newly  damaged  materials: 

- [eJ  mJ  + e”  (AV  Pj)l/Mj'^^ 


3.2  SYNCHRONIZATION  OF  MASS  TRANSPORT  AND  THE  MOVEMENT 

OF  DAMAGE  REGION  BOUNDARIES 

Previously  the  HELP  code  en^loyed  two  different  trans- 
port algorithms:  one  for  the  tracer  particles,  which  define 
the  material  interfaces  and  the  damage  region  boundaries  and 
another  for  the  mass.  As  described  in  Section  3.1.1  and 
illustrated  by  Figure  3.4a,  the  velocity  of  each  tracer 
particle  is  determined  by  an  interpolation  procedure  which 
area  weights  the  velocities  of  at  most  fout'  cells  in  the 
vicinity  of  the  particle,  and  which  thereby  define  a spatially 
continuous  velocity  field. 

The  original  HELP  algorithm  for  transport..ng  mass,  on 
the  other  hand,  presumed  a discontinuous  velocity  field,  i.e.. 
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a.  Tha  tracer  particle  velocity  is  based 
on  values  at  locations  denoted  by  x, 
which  were  interpolated  from  neigh- 
boring cell  centered  values. 


mm 

t 

b.  In  the  original  HELP  treat- 
ment, a single  transport 
velocity  (at  x)  computed  at 
each  cell  boundary. 


Figure  3.4. 


c.  In  the  modified  transport 
treatment,  transport  ve- 
locities Cnt  x)  are  com- 
puted at  three  points  for 
each  cell  boxmdary. 

velocity  Interpolations  used 


A comparison  of  the 
in  moving  the  tracer  particles  and  in  computing 
the  transport  quantities. 

‘ (1  • 


i 


28 


it  was  based  on  the  assiaroptlon  that  all  of  the  mass  in  each 
cell  had  the  velocity  of  the  cell  ceriter.  As  illustrated  by 
Figure  3.4b,  the  volume  of  material  transported  from  one  cell 
to  another  was  based  only  on  the  cell-centered  velocities  of 
the  two  cells. 


3.2.1  Use  of  Continuous  Velocity  Field  in  Mass  Transport 

To  synchronize  the  movement  of  the  tracer  particles 
and  the  mass,  the  calculation  of  the  transport  volumes  for 
transporting  mass  was  rewritten  so  as  to  be  based  on  a 
spatially  continuous  velocity  field.  Consider  the  mass  trans- 
port in  the  axial  direction  as  illustrated  in  Figure  3.4c. 

The  "map  line"  is  the  locus  of  points  such  that  each  point 
exactly  reaches  the  cell  interface  when  transported  through  the 
time  step  At  with  its  locally  interpolated  axial  velocity  com- 
ponent. If  every  material  point  in  the  grid  (together  with  its 
associated  momentum  and  energy)  is  moved  through  the  time  step 
At  with  its  local  (Interpolated)  axial  velocity,  and  if  each 
tracer  point  is  moved  likewise,  then  axial  mass  transport  is 
consistent  with  axial  tracer  motion.  The  same  argument  applies 
to  the  radial  transport  and  tracer  motion. 

If  a material  interface  intersects  a cell  boundary,  the 
transport  volume  determined  by  the  map  line  is  apportioned 
among  the  constituents  on  the  basis  of  the  fractional  areas 
subtended  by  each  material. 


3.2.2  Division  of  Transport  Phase  Into  Axial  and  Radial 

Subphas^ 

To  simplify  the  logic  of  the  new  transport  scheme, 
the  transport  phase  of  the  computational  cycle  has  been 
divided  into  two  subphases  with  only  axial  (or  radial)  trans- 
port taking  place  in  the  first  subphase  and  radial  (or  axial) 
transport  in  the  second.  To  prevent  bias,  the  order,  radial- 
axial  or  axial-radial,  is  reversed  each  computational  cycle. 
The  tracer  particles  as  well  as  the  mass  are  moved  in  these 
two  subphases. 


3 . 3 REVISED  COMPUTATION  OF  PRESSURES  AND  DEVIATCRIC 
STRESSES 

The  computation  of  the  deviatoric  stresses  and  the  cal- 
culation of  the  effects  of  both  the  pressure  and  the  deviatoric 
stresses  has  been  modified  in  the  NAG-FRAG  version  of  the  HELP 
code  as  described  in  the  sections  that  follow. 


3.3.1  Computinq  Separate  Stress  States  for  Materials 

IRHixerCeiTi ' 

In  order  to  apply  the  NAG-FRAG  failure  model  to  each 
material  in  mixed  or  free  surface  cells,  an  accurate  account 
of  the  stress  state  of  each  material  is  necessary.  Previously, 
tensile  pressures  in  all  interface  cells  2md  the  deviatoric 
stresses  of  free  surfaces  cells  were  set  to  zero,  thus  making 
it  impossible  to  apply  a failure  criterion  to  these  cells. 
Table  3.1  summarizes  the  changes  made  to  the  calculation  of 
the  pressures  and  deviatoric  stresses  in  mixed  and  free  sur- 
face cells  so  that  the  NAG-FRAG  failure  model  could  be  applied 
individually  to  the  materials  in  these  cells. 


3 V 3 . 2 Combining  the  Hydrodynamic  and  Strength  Phases 

of  HELP 

Previously  the  HELP  computational  cycle  was  divided 
into  three  phases,  the  hydrodynamic  phase  (HPHASE)  the  strength 
phase  (SPHASE) , and  the  transport  phase  (TPHASE)  . In  the 
hydrodynamic  phase  the  cell  pressures  were  computed  and  the 
effects  of  these  pressures  on  the  cell  velocities  and  internal 
energies  were  determined;  an  intermediate  set  of  cell  velocities 
and  energies  were  stored.  In  the  strength  phase  the  strain 
increments  were  computed  (using  the  post-HPHASE  intermediate 
velocities)  and  the  deviatoric  stresses  of  the  cells  were  up- 
dated accordingly.  These  updated  stress  deviators  were  then 
applied  in  the  momentum  and  energy  equations  and  another 
intermediate  set  of  cell  velocities  and  internal  energies 
were  defined. 

The  first  two  phases,  HPHASE  and  SPHASE  have  been  com- 
bined to  cor  form  with  the  requirements  of  the  BFRACT  failure 
criterion,  uhich  considers  the  total  stress  components  normal 
to  the  bln  -ingles,  and  to  facilitate  the  mapping  of  the  damage 
region  stress  states  into  the  Euleri2m  grid.  Thu..,  the  effects 
of  the  preisure  and  deviatoric  stresses  are  now  computed  in 
the  same  phase  (SPHASE  for  stress  phase) , and  only  one  set  of 
intermediate  velocities  and  Internal  energies  are  generated 
prior  to  the  transport  phase  (still  TPHASE)  . 


3.3.3  Mapping  the  Damage  Region  Stresses  into  the 
Eulerlan  Grid 


On  each  computational  cycle  the  Eulerlan  cell  stresses 
are  defined  by  the  following  procedure.  First  the  stresses 
of  each  damage  region  are  computed  by  the  NAG-FRAG  subroutines 
DFRACT  and/or  BFRACT.  Next  the  HELP  subroutines  determine  the 
stress  state  of  the  material  in  all  of  the  Eulerlan  cells. 
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Next,  if  the  stresc  state  of  undamaged  material  in  an  Eulerian 
cell  satisfies  the  NAG-FRAG  failure  condition,  a new  damage 
region  is  created  to  redefine  the  stress  state  of  this  newly 
formed  damage  region.  Finally,  the  stresses  c£  all  the  damage 
regions  (old  and  new)  are  mapped  into  the  Eulerian  grid,  thus 
altering  the  stress  state  of  those  cells  containing  damaged 
material.  For  exanqple,  the  final  cell  pressure,  P^,  of  a 
single  material,  pure  cell  is  defined  as  follows: 


if  the  cell  contains  only 
undamaged  material 


Z V* 


n 


if  the  cell  contains  only 
damaged  material 

r 


V*P®  + ^V^P^  if  the  cell  contains  both 

ry damaged  and  undamaged 

^ ^cell  material 


where 

V*  • the  volume  common  to  the  n damage  region  and  to 
the  cell 

V®  » the  volume  of  undamaged  material  in  the  cell 

£ ^ V 

P^  > the  pressure  in  the  n damage  region 

A 

P « the  pressure  of  the  undamaged  material  in  the 
cell. 

I 

The  volume,  V , coranion  to  a given  damage  region  and  a given 
cell  is  conduced  using  the  current  positions  of  the  d2unage 
region  tracer  particles  (which  define  the  boundary  of  the 
region)  and  the  fixed  boundaries  of  the  cell.  These  partial 
volumes  are  computed  every  computational  cycle  for  every 
cell  which  contains  the  boundary  of  one  or  more  damage  regions. 


4.  AUXILIARY  IMPROVEMENTS  TO  THE  HELP  CODE 


Three  auxiliary  improvements  were  made  to  the  HELP  code 
as  the  NAG-FRAG  failure  models  were  incorporated  and  tested. 
Options  were  added  for  making  the  top  and  right  grid  bo'indaries 
either  transmittive  or  reflective,  and  for  adding  either,  or 
both,  a linear  and  a quadratic  artificial  viscosity  term  to 
cell  boundary  stresses.  Further,  the  a-e  phase  transition  in 
the  Hugoniot  of  iron  in  compression  was  ii^plemented  into  the 
iron  equation  of  state.  These  improvements  will  be  described 
in  the  sections  that  follow. 


4.1  ADDITIONAL  REFLECTIVE  GRID  BOUNDARY  OPTIONS 

In  the  1975  documented  version  of  HELP,  the  left  grid 
boundary,  V7hlch  is  an  axis  of  symmetry  when  cylindrical  coordi- 
nates are  used,  is  always  reflective,  the  bottom  grid  boundary 
can  be  either  transmittive  or  reflective,  and  the  right  and 
top  grid  boundaries  are  always  transmittive.  (The  transmittive 
and  reflective  grid  boundary  conditions  are  described  in  detail 
in  Chapter  III  of  the  1975  documentation  of  HELP  [Reference  2].) 
The  tapered  flyer  plate  plane-strain  calculations,  which  were 
performed  in  order  to  test  out  the  implementation  of  the  NAG- 
FRAG  failure  models,  required  that  the  right  as  well  as  the 
left  grid  boundary  be  reflective.  Thus  the  option  to  have  the 
right  and  top  grid  boundaries  be  reflective  was  iiqplemented, 
even  though  the  reflective  boundary  condition  was  required  only 
on  the  right. 

In  the  new  version  of  KELP,  therefore,  the  user  dictates 
the  condition  of  the  bottom,  right  and  top  grid  bo'indaries  when 
defining  the  Z-block  variables  CVIS,  CVISR,  and  CVIST,  respec- 
tively, in  the  input  deck.  A zero  value  (which  is  the  default 
value)  indicates  a transmittive  condition;  whereas,  -1.  indicates 
a reflective  condition.  (Note:  The  opposite  convention  is 
used  in  the  other  versions  of  HELP,  wherein  -1.  indicates  a 
transmittive  boundary  condition.  Also,  the  location  of  CVIS 
in  the  Z-block  has  been  changed  from  27  to  125.) 


4.2  LINEAR  AND  QUADRATIC  ARTIFICIAL  VISCOSITY 

The  form  of  the  artificial  viscosity  added  to  cell 
boundary  stresses  in  the  NAG-FRAG  version  of  HELP  is  more  accu- 
rate than  that  used  in  the  standard  1975  version  of  the  code. 

The  expression  for  linear  artificial  viscosity  when 
applied,  for  example,  to  the  right  cell  boundary  is  of  the  form 
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- QLIN*c*p*/X*Kjy^ 

where 

QLIN  * linear^ artificial  viscosity  coefficient 
c ■ sound  speed 

r 

p « material  density 
A « cell  area  - AR*AZ 

■ normal  strain  rate  in  R direction. 


In  the  1975  documented  version  of  HELP,  this  expression  was 
simplified  by  assuming  AR  ■ AZ,  ■ 1 and  c « /P7f>  which 

led  to 

Q-  - /T"  * tu(K+l)  - u(K)]  . , 

L p 

I 

Such  a sound  speed  calculation  is  suitable  for  an  ibeal  gas 
(in  which  c • /y  P/p)  but  a poor  approximation  for  materials 
represented  by  a Tillotson  equation  of  state.  ' 

j 

The  linear  artificial  viscosity  used  in  the  NAG-FRAG 
version  (1)  does  not  require  square  zoning,  (2)  applies  a 
user  defined  coefficient  QLIN  * Z(93)  and  (3)  uses  a more 
accurately  calculated  sound  speed:  ' ' 

c^  for  a Tillotson  equation  of  state  (materials  1-19) 

for  an  ideal  gas  or  JWL  equation  of  state  i. 

^ (materials  20-25) 


In  addition,  a quadratic  artificial  viscosity  term, 
Q.,  may  be  added  in  compressing  regions: 

^ { 

Qq  ■ QQUAD*p*A*(ej^p)^ 


where  QQUAD  « Z(94)  is  defined  by  the  user.  (NOTE:  since 
the  use  of  artificial  viscosity  is  now  governed  by  the  input 
variables  QLIN  and  QQUAD,  the  formerly  used  input  flag  LVISC 
has  been  deleted  from  the  Z-array, ) 


4.3  IMPLEMENTATION  OF  THE  a-C  PHASE  CHANGE  IN  IRON 

In  addition  to  the  specific  work  involving  the  frag 
mentation  models,  the  equation  of  state  for  iron  has  been 


modified  to  account  for  the  a-e  phase  change  that  occurs  at 
pressures  in  the  vicinity  of  130  kbars.  This  modification 
was  made  since  any  failure  model  depends  on  dynamic  quantities 
calculated  from  the  constitutive  relation,  and  the  accuracy 
with  which  BFRACT  and  DFRACT  predict  fragmentation  in  rolled 
homogeneous  armor -and  penetrator  steels  is  therefore  highly 
dependent  upon  the  accuracy  of  the  constitutive  model  over  the 
entire  range  of  stresses  and  pressures  expected. 

The  usual  form  for  the  equation  .of  state  employed  in 
the  HELP  code  is  that  due  to  Tillotson^^^  modified  to  give  a 
smooth  transition  between  condensed  and  expanded  states.  For 
t>e  condensed  states,  i.e.,  when  p/p.  > 1,  or  for  any  cold 
states,  E < E the  equation  of  state*^has  the  form 
® - b) 


P - P, 


a + 


+ 1 


Ep  + Ap  + Bp' 


(1) 


For  expanded  hot  states,  i.e.,  when  p/p  < 1 and  E < E^,  the 
equation  of  state  has  the*  form  ' ! 


p « p - aEp  + 


•¥  1 


+ Ape 


-e(v/v^-l) 


-a(v/v  -1) 

e ° . (2) 


A smooth  transition  between  the  condensed  and  expanded  states 
is  insured  by  a transition  equation  for  the  intermediate 
region  defined  by  E < E < E'  and  p/p_  < ‘1.  This  blended 

portion  of  the  equation  of  state  has  the  form 

1 , «■ 


(E  - E 


s^P 


E * 


- E)Pe 


'.'b' 


(3) 


1 

In  Equations  (1)  through  (3)  p,  E and  p are, pressure,  specific 
Internal  energy,  and  mass  density,  respectively,  n ■ p/Pq  " 
p + 1 « V /v,  and  p , a,  b,  E , E , E',  A,  B,  o^and  0 are  con- 
stants for  the  particular  material.  , 


The  material  parameters  employed  in, the  Tlllotson 
equation  of  state  for  iron  arc  given  in  Table  4.1. 


The  Tillotson  equation  of  state  for  compressed  states 
of  iron  may  also  be  written  as 


TABLE  4.1.  MATERIAL  PARAMETERS  FOR  IRON 


, ^><^TERIAL 

property"^ 


1 pQ  (g/cm^) 

v7.856  . 

2 

A (dyntts/cm  ) 

Sea  Text 

2 

B (dyn«s/cn  ) 

Saa  Text 

a 

0.5 

b 

1.5 

1 

a 

T * r 

5.0 

B 

5.0 

Eq  (ergi/g) 

9.5  X 10^® 

E^  (arga/g) 

2.44  X 10^® 

EJ  («rg«/g) 

lw02  X 10^^ 

2 

G (dynes/cm  ) 

f.O  X 10^^ 

Yq  (dynari/cm^) 

6.0  X 10^ 

(<^lyne  1/051^) 

0 

Yj  (dyn«s/cm^) 

0 

(«rgf/g) 

1.3  X 10^® 

1 

2 

S*  (dyn««/oin  ) 

3.8  X 10^° 

** 

1 .• 

Spall  atrangth  > 

\ 

• 

'■’i 
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(4) 


for  n ^ 1*0  is 


Ejj  . / (5) 


The  second  term  in  the  expression  for  f(p)  insures  that  the 
pressure  P will  be  equal  to  the  Hugoniot  pressure « P.,«  when 
the  energy  E is  equal  to  the  Hugoniot  energy,  'E^^. 

The  analytic  fit  to  the  Hugoniot  of  iron  in  con5)ression 
as  given  by  Mader^^^  was  modified  for  use  in  the  HELP  code. 

The  a-e  phase  transition  is  bridged  by  a two-piece  analytic 
fit  which  mairtains  continuity  of  the  pressure  and  its  first 
derivative  along  the  Hugoniot  curve.  The  resulting  approxi- 
mate equation  of  state  assumes  that  the  phase  change  is 
instantaneous. 

For  low  pressure  states  (0  < Pu  < 125  kbar)  the  analytic 
fit  to  the  Hugoniot  is 


5 , (6) 

(1  - 

where 

3 

« 7.856  gm/cm 
■ 4.63  X 10^  cm/sec 
« 1.275 

and 

e * 1 - Pq/p  . 


P - f(p)  + 


a + 


1 + 


5 


pE 


where  the  appropriate  form  of  f(p). 


f{p)  - Ph  ■ 


a + 


1 + 


'H 
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Equation  (6)  for  the  pressure  In  dyne/cm  is  valid  for  the 
range  of  specific  volumes 


0.11929584  < V < 0.12729124  - l/p^ 

For  the  pressure  range  from  125  kbar  to  140  kbar,  the 
Hugonlot  pressure  is  fitted  as 

Pjj  » A + fl  (AV)  + (AV)^  ♦ (AV)^®  (7) 

( 

j 

where  > 

A - 140  X 10^ 

B - - 4.0  X 10^^ 

- 2.383828348  x 10^^ 

- - 6.060826938  x 10^^ 

AV  « V - 0.1162 

I I -1 

and  the  range  of  specific  volumes  is 
0.1162  < V < 0.11929584 

lur  the  pressure  range  from  140  kbar  to  155  kbar,  the 
Hugonlot  pressure  is  given  by 

Py  - A + B (AV)  + (AV)^  + (AV)^°  (8) 

with  A and  B as  before  and 

Cj^  - - 7.090650579  x 10^^ 

- 1.224016343  x 10^^ 

AV  £ V - 0.1162 

The  range  of  specific  volomes  is 


0.11335815  < V < 0.1162 


For  the  pressure  range  155  to  364  kbar,  the  analytical 
fit  to  the  Hugoniot  is 


with 

3 

Pq  » 7.856  gn/cm 

« 3.20  X 10^  cm/sec 
Sjj  - 2.25 
e ■ 1 - Pq/p  , 

and  the  range  of  specific  volumes  is 


(9) 


0.1055  < V < 0.11335815 

( 

For  pressures  from  364  kbar  to  457  kbar  the  following 
equations  are  suggested. 

For  volumes  in  the  range 
0.1030  < V < 0.1055 

2 

the  Hugoniot  pressure  in  (dyne/cm  ) is 

2 3 

Pjj  - A + B .AV)  + C (AV)^  + D (AV)-*  (10) 

with 

A - 3.6432733  x 10^^ 

B - - 3.7667780  x 10^^ 

C - - 1.712352  x 10  ‘ 

D - 3.6615110  X 10^® 

AV  ■ V - 0.1055 
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For  specific  volumes 

r .X  ‘‘ 

V < 0.1030 


the  Hugoniot  pressure  is 


u 


- Po^o 


(1  - Se) 


(11) 


with 


■ 7.856  gm/cm 

■ 3.78  X 10^  cm/sec 

S - 1.652 
e ■ 1 - Pq/p 


•J  f 


For  presHures  in  the  range  0 ^ P < 125  kbar,  or  equivalently 
for  volumes  in  the  range 


0.11929584  ^ V 0.12729124 


•.  ) 


r t j i<  * 


the  Hugoniot  energy  is  computed  with 


? j 


=H  - I '■h  <''o  - V)  . 


(12) 


where 


- 0.12729124 
o 


r ; 


For  the  volume  range 


0.1066867  < V < 0.11929584 


the  energy  is  computed  from 


- J '■l  '''o  - ''l>  * ? •'■h  ♦ ''l>  'Vl  - V) 


(13) 


* *-  ^ * 
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with 


1.25  X 10^^  dyne/cni* 


/ 


V,  • 0.11929584  cnVgm 


/ 


For  the  volume  lees  than  0.1066867,  Equation  (12)  is 
used  to  compute  the  Hugoniot  energy.  / ! 


/ 


; J- 


The  modified  Hugoniot  for  modeling  the  phase  change  in 
iron  is  plotted  in  Figure  4.1  for  the  pressure  range  between  0 
and  400  kbar  and  a magnified  plot  of  tM  fits  to  the  transition 
region  is  given  in  Figure  4.2.  i 


For  condensed  iron  at  pressures  less  than  0 kbar  the 
form  of  the  f(p)  equation  is  i 


/ 


f{p)  ■ Au  + Bu‘ 


(14) 


where 


A 

B 


1.28  X 10 


12 


1.05  X 10 


12 
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Figure  4.1,  Hugoniot  for  iron. 
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5.  APPLICATIONS  OF  THE  NAG-FRAG  DUCTILE 


AND  BRITTLE  FAILURE  MODELS 


The  NAG-FRAG  version  of  the  HELP  code  has  been  applied 
to  three  situations: 

1.  an  1145  aluminum  tapered  flyer  plate  impacting 
a target  of  the  same  material, 

2.  an  iron  tapered  flyer  plate  impacting  a target 
of  the  same  material,  and 

3.  the  fragmentation  of  a cylindrical  iron  shell 
containing  a detonated  high  explosive  (COMP-B) . 

The  aluminum  was  treated  as  a ductile  material  in  the 
first  calculation  and  the  iron  as  a brittle  material  in  the 
second  and  third.  Thus  both  of  the  NAG-FRAG  models  BFRACT  and 
DFRACT  were  exercised. 

The  results  obtained  from  these  calculations  were  then 
compared  with  those  obtained  by  SRI  with  the  Lagrangian  HEMP 
and  PUFF  codes. 


5.1  1145  ALUMINUM  TAPERED  FLYER  PLATE  IMPACTING  AN 

1145  ALUMINUM  TARGET 

A two-dimensional  Lagrangian  calculation  simulating  the 
impact  of  an  1145  aluminum  tapered  flyer  plate  on  a target  of 
the  same  material  usina  the  ductile  failure  model  DFRACT  has 
been  reported  by  SRI.‘^'  This  calculation  has  been  repeated 
using  the  same  failure  model  in  the  HELP  code. 

i 

A description  of  the  target  and  flyer  plate  configuration 
along  with  the  zoning  used  in  both  the  Lagrangian  and  HEMP 
calculations  are  shown  in  Figure  5.1.  The  material  properties 
and  fracture  parameters  used  in  the  calculation  are  given  in 
Table  5.1. 

In  the  HELP  calculation  several  options  were  considered 
concerning  the  treatment  of  the  target-flyer  plate  interface: 

(1)  not  allowing  tensile  pressures  in  interface  cells,  (2)  allow 
ing  tensile  pressures  but  not  allowing  separation  at  the  inter- 
face and  (3)  removing  the  flyer  plate  when  the  interface  first 
becomes  tenjile.  It  was  found  that  the  third  option  gave  the 
closest  agreement  with  the  Lagrangian  calculation.  (In  the 
Lagrangian  calculation  interface  separation  was  permitted, 
hence  this  option  appears  the  most  appropriate.)  Also,  several 
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TABLE  5.1.  MATERIAL  CONSTANTS  AND  FRACTURE 
PARAMETERS  FOR  1145  ALUMINUM 


Til : 


-,14  r-  ■ 
; ! ■ 


SYMBOL 


VALUE 


DEFINITION 


760  kbart 


1500  kbars 


2.04 


. . i . ■ ■ t ! 

Coefficients  in  the  Mie-Gruneisen 
equation  of  state  , ! 

■ M 

p « (Am’ + bm^  Vcm^XI  - j w)  + rpE 

I ■ ' 

M ■ <P/Po  ‘ 


300  kbars 


Shear  modulus 


1.1547  kbars 


2.70  gm/cm' 


Yield  strength  in  shear 
Initial  density 


-0.01  cm  /dyne/sec 


Void  growth  constant 


-4.0  kbars 


Void  growth  threshold 


lO”^  cm 


Void  nucleation  radius  parameter 


9 3 

3.0  X 10  /cm  /sec 


-3.0  kbars 


-0.4  kbars 


Parameters  for  void  nucleation 


N - "o  * 


•'■-"no'/'’! 


■%T,  ^ ^ 


levels  of  linear  artificial  viscosity  were  tested  in  the  HELP 
code.  Although  the  results  were  generally  smoother,  the  damage 
levels  were  significantly  reduced  when  the  artificial  viscosity 
was  added.  Accordingly,  in  the  HELP  calculation  whose  results 
are  presented  here  no  artificial  viscosity  was  added  and  the 
flyer  plate  was  removed  at  about  0.71  ps. 

Since  the  rate  of  both  void  nucleation  and  growth  is  a 
strong  function  of  tensile  stress,  any  spurious  behavior  in 
the  stress  field  might  well  lead  to  a prediction  of  very  high 
localized  damage.  On  the  other  hand,  while  artificial  viscosity 
is  effective  in  damping  out  such  behavior,  it  was  found  in  the 
test  case  to  signif iciantly  alter  the  damage  throughout  the 
entire  grid.  Experience  suggests  a reasonable  means  of  treating 
this  problem  might  be  to  require  the  tensile  stress  to  exceed 
the  nucleation  threshold  for  several  cycles  before  allowing 
material  to  fail. 

A measure  of  the  dynamic  response  of  the  target  free 
surface  to  flyer  plate  impact  is  given  by  Figure  5.2  in  which 
the  average  velocity  across  this  surface  is  plotted  as  a 
function  of  time  for  both  the  Lagrangian  and  HELP  calculations. 

Two  measures  of  target  damage  level,  the  void  volume 
fraction  and  the  number  of  voids  per  cm^,  as  computed  by  the 
HELP  code  are  shown  in  Figures  5.3.  and  5.4  at  time':i  0.93,  1.17, 
and  1.62  ys.  (Corresponding  values  of  these  quantities  ob- 
tained from  the  HEMP  calculation  at  1.46  ms  are  shown  in 
Figures  5.5  and  5.6).  It  may  be  noted  that  the  areal  extent 
of  damage  increases  little  after  0.93  ms.  However,  the  level 
of  d2unage  increases  due  to  both  nucleation  and  growth  of 
voids.  The  damage  levels  existing  in  both  the  SRI  and 
calculations  after  target  recompression,  during  which  the 
damage  level  remains  constant,  is  shown  in  the  surface  plots 
of  Figures  5.7  and  5.8. 


5.2  ARMCO  IRON  TAPERED  FLYER  PLATE  IMPACTING  ON  ARMCO 

IRON  TARGET 

SRI  reported  a second  2-D  Lagrangian  impact  calcu- 
lation simulating  the  impact  of  an  Armco  iron  tapered  flyer 
plate  on  a target  of  the  same  material  using  the  brittle 
failure  model  BFRACT.  This  calculation  was  repeated  using 
the  same  failure  model  in  the  HELP  code. 

A description  of  the  target  and  flyer  plate  configura- 
tion along  with  the  zoning  used  in  both  the  Lagrangian  (HEMP) 
and  Eulerian  (HELP)  calculations  are  shown  in  Figure  5.9.  The 
material  properties  and  fracture  parameters  used  in  the  cal- 
culation appear  in  Table  5.2. 
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Figure  5.3.  Relative  void  “olume  of  damage  regions  at  three 
times  after  impact  in  the  1145  aluminum  flyer 
plate  calculation.  (Numbers  shown  are  to  be 
multiplied  by  10”‘,>  Impact  plane  is  at  bottom, 
and  thick  end  of  flyer  plate  is  at  right. 
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Figure  5.4. 


Humber  of  voids  per  cubic  centimeter  in  damage 
regions  at  three  times  after  impact  in  the  1145 
aluminum  flyer  plate  calculation.  (Humbcrs 
shown  are  to  be  multiplied  by  10^.)  Impact 
plane  is  at  bottom,  and  thick  end  of  flyer  plate 
is  at  right. 
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Figure  5.5  Relative  void  volume  of  material  in  the  1145 
aluminum  target  at  1.46  ys  after  Impact  as 
computed  by  the  Lagrangian  (HEMP)  code,  (Num- 
bers shown  are  to  be  multiplied  by  10'°.). 
Impact  plane  is  at  bottom,  and  thick  end  of 
flyer  plate  is  at  right. 
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Figure  5.6  Number  of  voids  per  cxibic  centimeter  of  material 
in  the  1145  aluminum  target  at  1.46  ys  after  im- 
pact as  computed  by  the  Langrangian  (HEMP)  code. 
(Numbers  shown  are  to  be  multiplied  by  10^.). 
Impact  plane  is  at  bottom,  and  thick  end  of  flyer 
plate  is  at  right. 
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Eulerian  (HELP)  Calculation 


Figure  5.8.  Nuaber  of  voids  per  c\ibic  centimeter  in  the  1145  aluminun  target 
after  reconpression  as  computed  by  both  the  Lagrangian  (HEMP) 
code  and  the  Eulerian  (HELP)  code. 
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TABLE  5.2.  MATERIAL  CONSTANTS  AND 
FRACTURE  PARAMETERS  FOR  ARMCO  IRON 


SYMBOL 


VALUE 


DEFINITION 


A 

B 

C 

r 

G 

Y 


3 


go 

% 

°no 

TSR(IC) 

TSR(ll) 

TSRa2) 

TSR(13) 


1589  kbars 
5170  kbars 
51700  kbars 
1.69 

819  kbars 
6.062  kbars 
7.85  gm/cm^ 

2 

-0.00012  (r.i^  /dyne/sec) 
0.0 
0.004 


''  Coefficients  in  the  Mie- 
Gruneisen  equation  oi  state 

P 


(Am  + bm^  cm^)(1  - j )+rpE 


5.0  X 10 


8 


-11.2  kbars 

-0.74  kbars 

0.25 

1.0 

0.2 

1.0 


(p/Pq  - 1) 


Shear  modulus 
Shear  yield  strength 

Initial  density 

1 

Void  growth  constant 

Void  growth  threshold 

Void  nucleation  radius 
parameter 

Parameters  for  void  nucleation 


y N « e 


)/a, 

(Pil/  no  1 


Fragmentation  parameters 
provided  by  SRI 


As  in  the  aluminum  flyer  plate  calculations,  the  iron 
flyer  plate  was  removed  when  the  projectile-target  interface 
first  became  tensile  in  both  the  HEMP  and  HELP  calculations. 

The  separation  occurred  at  about  1.32  us. 

To  get  the  best  agreement  between  the  HEMP  and  HELP 
results,  both  linear  and  quadratic  artificial  viscosity  terms 
were  added  to  the  ceil  boundary  stresses  in  compressing  regions 
of  the  HELP  grid;  in  the  expanding  regions  only  the  linear 
viscosity  term  was  added.  The  form  of  the  artificial  viscosity 
in  this  version  of  HELP  is  described  in  Section  4.2.  The  linear 
and  quadratic  coefficients  used  in  this  calculation  were  0.1 
and  4.0,  respectively. 

Three  measures  of  damage  in  the  brittle  target  material, 
the  number  of  cracks  per  cm^,  the  void  volume  fraction,  and 
the  fraction  of  fragmented  material,  as  computed  by  HELP,  are 
shown  in  Figures  5.10,  5.11  and  5.12  at  times  1.18,  1.46,  and 
2.11  us.  (Corresponding  values  of  the  number  of  cracks  per 
cm3  and  the  fraction  of  fragmented  material  as  computed  by  the 
HEMP  code  are  shown  In  Figures  5.13  and  5.14  at  1.24,  1.51  and 
2.33  U8«  The  void  volume  fraction  was  not  an  edited  quantity 
in  the  HEMP  calculations.)  It  is  evident  from  these  figures 
that  failure  begins  at  the  end  of  the  plate  impacted  by  the 
thin  part  of  the  flyer.  This  is  expected  (and  actually  occurs 
in  the  aluminum  impact  calculation  as  well)  since  the  rarefac- 
tion wave  from  the  free  surface  of  the  flyer  has  less  far  to 
travel  at  the  thin  end.  Although  the  damage  first  appears  at 
the  so-called  thin  ond,  the  greatest  damage  is  ultimately  in- 
curred at  the  end  impacted  by  the  thick  part  of  the  flyer. 

These  figures  also  reveal  a significant  increase  in  the  damage 
and  fragmentation  of  the  target  near  the  thick  end  of  the  flyer 
between  1.46  and  2.11  us;  whereas,  the  damage  at  the  Lhin  end, 
due  to  both  nucleation  and  growth  of  cracks,  does  not  increase 
significantly  after  1.46  us. 

The  surface  plots  of  Figures  5.15  and  f>.16  show  the 
nu:nber  of  cracks  and  the  amount  of  fragmentation  predicted 
by  the  Lagrangian  (HEMP)  code  and  the  Eulerien  (HELP)  code  at 
2.33  and  2.11  us,  respectively.  Although  the  HELP  calculation 
predicts  a slightly  larger  number  of  cracks  per  cm^  near  the 
thin  end  of  the  flyer,  and  somewhat  more  damage  near  the  back 
free  surface  of  the  target,  the  regions  of  complete  fragmentation 
and  the  concentration  of  cracks  at  the  high  damage  end,  as 
predicted  by  the  two  codes,  are  in  very  good  agreement. 


5.3  FRAGMENTING  MUNITIONS 

A 1-D  calculation  was  performed  with  the  NAG-FRAG 
version  of  HELT‘  in  which  the  expansion  and  subsequent  failure 
of  an  infinitely  long  iron  cylinder  filled  with  detonated  COMP  B 
explosive  was  simulated.  In  order  to  verify  the  HELP  prediction 
of  failure  near  the  free  surface  of  the  cylinder,  the  results 
were  compared  with  those  obtained  by  SRI  using  the  1-D  PUFF  code. 
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Figure  5.10.  Number  of  cracks  per  cubic  centimeter  in  damage 
regions  at  three  times  after  impact  in  the  Armco 
iron  flyer  plate  calculation.  Impact  plane  is 
at  bottom  and  thick  end  of  flyer  is  at  left. 
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Figure  5.11. 


Void  volume  fraction  in  damage  regions  at  three 
times  after  impact  in  the  Armco  iron  flyer  plate 
calculation.  Impact  plane  is  at  bottom  and  thick 
end  of  flyer  is  at  left.  Values  of  -2.0  correspond 
to  regions  that  are  completely  fragmented. 
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Figure  5.12 


. The  fraction  of  fragmented  material  in  damage 
regions  at  three  times  after  impact  in  the 
Armco  iron  flyer  plate  calculation.  Impact 
plane  is  at  bottom,  and  thick  end  of  flyer  is 
at  left. 
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Figure  5.13  Number  of  cracks  per  cubic  centimeter  of  material 
in  the  Armco  iron  target  at  three  times  after  im- 
pact as  computed  by  the  Lagrangian  (HEMP)  code. 
Impact  plane  is  at  bottom,  and  thick  end  of  flyer 
plate  is  at  left. 
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Figure  S.14  The  fraction  of  fragmented  material  in  the  Armco 
iron  target  at  three  times  after  impact  as  com- 
puted by  the  Lagrangian  (HEMP)  code.  Impact  plane 
is  at  bottom  and  thick  end  of  flyer  plate  is  at 
left. 
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Lagrangian  (HEMP)  Calculation  Eulerian  (HELP)  Calculation 
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Lagrangian  (HEMP)  Calculation  Eulerian  (HELP)  Calculation 
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The  dlmenslonB  of  the  explosive  and  iron  and  the  grid 
configurations  used  in  the  PUFF  and  HELP  calculations  appear 
in  Figure  5.17  (HELP  requires  a ninirouro  of  three  rows  in 
the  grid.)  The  suiterial  constants  for  the  iron  are  given  in 
Table  5.3.  In  both  calculations  the  COMP  B explosive  was 
detonated  simultaneously  at  time  t > 0 by  setting  the  initial 
specific  internal  energy  of  all  explosive  cells  to  4.469  x 10^^ 
ergs/g.  The  pressure  of  the  COMP  B was  derived  from  a gamma- 
law  equation  of  state,  with  y **  2.841. 

''  The  iron  begins  to  fail,  in  both  calculations,  at  the 
free  surface  at  about  6 ps.  Both  codes  predict  significant 
damage  to  the  outer  ring  of  the  cylinder.  At  9.6  ms  HELP 
predicts  that  the  damage  extends  into  the  cylinder  wall  0.745  cm 
from  the  free  surface;  at  11.7  ys  PUFF  predicts  the  damaged 
material  extends  0.826  cm  into  the  cylinder  from  the  free  sur- 
face. 


Figures  5.18  through  5.20  show  the  number  of  cracks, 
the  average  crack  length  and  the  fraction  of  fragmented  materia] 
in  each  damage  region  as  computed  by  HELP  at  three  times  after 
the  cylinder  begins  to  fail. 


Figure  5.17.  Dimensions  and  grid  configuration  for  the  fragmenting  round 
calculation. 


TABLE  5.3.  MATERIAL  CONSTANTS  AND  FRACTURE 
PARAMETERS  FOR  IRON  CYLINDER 


SYMBOL  VALUE  DEFINITION' 


A 

1589  kbars  . 

Coefficients  in  the  Mie- 
Gruneisen  equation  of  state 

B 

, 5170  kbar* 

5 3 r 

p » (Am  + bm  + Cm  ) (i-jM)+rpE 

C 

51700  >.bar> 

M * (p/pQ  “ 

r 

1.69 

G 

772.9  kbars 

Shear  modulus 

Y 

5.774  kbars 

Yield  strength  in  shear 

P 

o 

7.85  gm/cm^ 

Initial  density 

3 

-5.5  X 10“^ 

2 

cm  /dyne/sec  Void  growth  constant 

» 

go 

-0.10  kbars 

Void  growth  threshold 

«n 

4 X lO”^  cm 

Void  nucleation  radius  parameuer 

K 

8 3 

6.0  X 10  /cm  /sec 

i 

t 

-9.0  kbars 

Parameters  for  void  nucleation 

N . N e 

no 

-1.3  kbars 

o 

- - ^ — -i 

I 
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t = 9.6  ys 


Figure  5.20.  Fragmentation  fraction  as  computed  by  the  HELP 
code  at  three  different  times  after  detonation 
in  the  1-D  fragmenting  round  calculation. 
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APPENHIX  A 

USER'S  GUIDE  TO  THE  VERSION  OF  HELP 
WHICH  EMPLOYS  THE  NAG-FRAG  FAILURE  MODEL 


The  documentation  presented  in  this  Appendix  serves 
as  a supplement  to  the  1975  documentation  of  the  HELP  code 
(Reference  2) . Diagrams  representing  the  organizational  mod- 
ifications as  well  as  descriptions  of  new  subroutines  and 
new  input  variables  are  presented  in  the  sections  that  follow. 
Additionally,  this  writeup  provides  an  algorithm  for  esti- 
mating the  core  storage  requirements  of  a given  configuration. 


A.l  ORGANIZATION  MODIFICATIONS  TO  THE  HELP  CODE 

The  motivation  for  reorganizing  parts  of  the  HELP  code 
has  been  given  in  Section  3 of  this  report.  As  indicated 
there,  the  division  into  three  phases,  a hydrodynamic  phase, 
a strength  phase,  and  a transport  phase  has  been  reduced  to 
two  phases  by  combining  the  effects  of  t cell  pressures  and 
the  cell  deviator  stresses.  Also,  the  Eulerian  cell  stresses 
are  modified  to  reflect  the  existence  of  damaged  materials  by 
a procedure  which  maps  the  stresses  of  the  Lagrangian  damage 
regions  into  the  Eulerian  grid.  Finally,  the  transport  phase, 
including  the  motion  of  the  tracer  particles,  is  executed  in 
two  directional  passes.  These  modifications  can  be  seen  in 
Figure  A.l  by  comparing  the  flow  diagram  of  the  1975  documented 
HELP  code  with  the  flow  diagram  of  the  present  NAG-FRAG  version. 


A. 2 NEW  SUBROUTINES 

In  addition  to  the  BFRACT  and  DFRACT  subroutines  pro- 
vided by  SRI,  several  additional  subroutines  were  created  to 
accomplish  the  modifications  mentioned  above  as  well  as  to 
generate  and  update  Lagrangian  damage  regions.  These  added 
subroutines  are  listed  in  alphabetical  order  and  briefly 
described  in  the  paragraphs  below.  Detailed  flow  diagrams  of 
SCALC,  MIXDMR  and  PURDMR,  the  major  routines  which  process 
the  damage  regions,  appear  in  Figures  A. 2,  A. 3,  and  A. 4. 

ADJID  (called  from  MIXDMR) 

Subroutine  ADJlD  is  called  in  place  of  CALFRC,  VOLFND 
and  INTSAV  in  1-D  calculations  (JMAX=3)  to  enlarge  and  reduce 
existing  damage  regions. 

BFRACT  (called  from  DAMAGE) 

Subroutine  BFRACT  was  provided  by  SRI  to  compute  the 
pressure,  deviator jc  stresses  and  damage  parameters  of  brittle 
materials  that  have  satisfied  the  BFRACT  failure  condition. 
Chapter  2 of  this  report  describes  the  BFRACT  failure  con- 
ditions as  well  as  the  algorithms  for  crack  nucleation  and 
growth. 
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1975  DOCUMENTED  HELP  HELP  WITH  NAG-FRAG  FAILURE  MODEL 


DAMAGE  (called  from  MIXDMR,  PURDMR  and  SCALC) 


This  subroutine  is  called  from  SCALC  each  computational 
cycle  (provided  at  least  one  damage  region  has  been  generated) 
to  update  the  properties  of  all  damaged  material.  It  is  also 
called  from  MIXDMR  and  PURDMR  whenever  a new  damage  region 
is  created  or  an  existing  damage  region  is  enlarged. 

Subroutine  DAMAGE  computes  the  strain  increments,  a 
rotation  increment,  a new  relative  volume  (or  density)  and 
an  estimated  new  internal  energy  for  the  damage  region(s). 

These  quantities  are  used  by  the  NAG-FRAG  models  to  update 
the  pressure,  deviator ic  stresses,  number  and  size  of  voids 
or  cracks  and  fragmentation  data  which  are  then  mapped  into 
the  Eulerian  grid. 

DFRACT  (called  from  DAMAGE) 

Subroutine  DFRACT  was  provided  by  SRI  to  compute  the 
pressure,  deviatoric  stresses  and  damage  pareuneters  of 
ductile  materials  that  have  satisifed  the  DFRACT  failure 
condition.  Chapter  2 of  this  report  describes  the  DFRACT 
failure  condition  as  well  as  the  algorithms  for  void  nuclea- 
tion  ani  growth. 

DVOLID  (called  from  CALFRC) 

Subroutine  DVOLID  ir  called  instead  of  VOLFND  in  1-D 
calculations  (JMAX*3)  to  compute  the  volume  of  a cell  that 
belongs  to  a given  damage  region. 

DVSTRS  (called  from  SCALC) 

Subroutine  DVSTRS  updates  the  cell-centered  deviatoric 
stresses  of  the  pure  Eulerian  cells  and  the  deviatoric  stresses 
of  each  material  in  the  multimaterial  Eulerian  cells.  Basically, 
this  subroutine  replicates  the  first  half  of  the  SPHASE  sub- 
routine in  the  unmodified  HELP  code.  However,  unlike  the 
earlier  version  of  HELP,  the  deviatoric  stresses  of  each 
constituent  of  a mixed  cell  are  computed  separately;  (previously, 
the  yield  strengths  and  shear  moduli  of  the  constituents  were 
averaged,  and  the  same  average  stress  state  was  associated 
with  all  of  the  materials  in  a mixed  cell)  . A better  definition 
of  the  stresses  of  the  constituents  in  a mixed  cell  was 
necersitated  by  the  separate  application  of  the  NAG-FRAG 
failure  condition  tc  each  constituent  of  a mixed  cell. 
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EDOT  (called  from  DVSTRS) 


Subroutine  EDOT  computea  the  spatial  derivatives  of 
the  velocity  field  ^ and  used  In  the  computation 

of  the  devlatorlc  strain ‘’incrtinents  for  a given  pure  cell  or 
for  one  of  the  materials  In  a given  multimaterial  cell.  Near 
the  material  Interfaces,  the  difference  equations  are  modi'* 

! f led  so  that  the  derivatives  are  based  only  on  the  velocity 
field  within  the  given  'material  package,  as  Illustrated  below: 


• Pure  cells: 


KL 

K 

^KL 

^KR 

^KR~^KL 

^kr’^kl 


1 


> -I 


Interface  cells: 


Material  A 


Material  B 


Xk-=‘kl 


' iX'  1 

INTSAV  (called  from  VOLFND  and  CALFRC) ' ‘ ' ‘ 

Subroutine  INTSAV  Is  called  when^  an  existing  damage 
region  is  being  reduced  or  enlarged,  i.e.,  when  the  undamaged 
material  In  a partially  damaged  cell  has  satisfied  the  failure 
criterion: 


r;n  large 
r.xistinq 
Damage 
Region 


Reduce  Existing 
Damage  Region; 
Create  a Now 
Region 


'.'hen  a damage  region  is  being  either  enlarged  or  reduced, 
INTSAV  removes  the  damage  tracers  that  are  interior  to  the 
failed  cell  and  adds  the  appropriate  intercepts  and  corner 
points  which  have  been  defined  by  VOLFND  and  stored  in  the 
XT,  YT  arrays.  The  identification  of  the  beginning  and  end- 
ing points  (IBGN,IEND)  of  each  subsequent  damage  region  are 
adjusted  to  account  for  the  tracers  added  and  deleted  in  this 
process. 

The  intercepts  of  the  reduced  damage  region  with  the 
failed  cell  are  saved  by  subroutine  INTSAV  and  added  by  sub- 
routine NWTCR  to  the  string  of  tracers  defining  the  new  damage 
region. 

Finally,  when  a single  material  in  an  interface  cell 
fails,  subroutine  INTSAV  stores  - in  the  XTS,YTS  arrays  - 
the  material  interface  tracer  coordinates  and  bhe  included 
corner  points  as  defined  by  VOLFND  and  CALFRC. 

INTSAV  is  called  from  subroutine  CALFRC  only  if  the 
region  being  enlarged  lies  entirely  inside  the  failed  cell. 

In  this  case  CALFRC  defines  the  boundary  of  that  region  to 
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be  coincident  with  the  cell  boundary  by  placing  five  of  its 
tracers  at  the  cell's  corners.  INTSAV  removes  the  remaining 
tracers  for  that  region  beyond  the  five  defined  by  CALFRC 
and  adjusts  the  counters  of  the  subsequent  regions  accordingly. 

MIXDMR  (called  from  SCALC) 

,/ 

When  a single  material  in  an  interface  ("mixed")  cell 
satisfies  the  NAG-FRAG  failure  criterion,  subroutine  MIXDMR 
is  called  either  to  create  a new  damage  region  or  to  enlarge 
an  existing  region.  In  either  case  one  or  more  existing 
damage  regions  may  be  reduced. 


Enlarge  an 
Existing 
Damage 
Region 


Create  a Nc\/ 
Region;  Reduce 
an  existing 
Region 


The  enlargement  or  reduction  of  an  existing  damage 
region  is  accomplished  by  adding  and/or  deleting  appropriate 
damage  tracer  particles,  recomputing  the  region's  volume  and 
mass,  and, in  the  case  of  enlargement,  averaging  the  stresses 
and  failure  parameters  of  the  added  material  with  those  of 
the  region  before  it  was  enlarged. 

The  tracer  strings  of  the  modified  regions  are  rede- 
fined by  making  use  of  the  mechanisms  of  CALFRC  and  VOLFND  to 
identify  the  intercepts  of  the  damage  region  boundary  with  ‘ 
the  boundaries  of  the  cell  and  to  identify  which  cell  corners 
lie  in  the  interior  of  the  region.  Subroutines  CALFRC,  VOLFND 
and  INTSAV  examine  the  value  of  the  variable  NCALL  to  discern 
whether  the  N^"  region  is  being  reduced  or  enlarged. 


In  MIXDMR  a region  is  enlarged  to  encompass  only  a 
single  material  in  a cell,  i.e.,,  not  the  whole  cell.  The 
procedure  for  doing  so  is  broken  down  into  two  steps,  each 
requiring  a call  tp  CALFRC  and  VOLFND.  ^ First,  the  damage  i 

tracers  interior  to  the  cell  are  removed  and  the  cell  corner 
points  and  intercepts,  identified  by  CALFRC  and  VOLFND,  are  /') 
added  by  INTSAV  as  if  the  entire  cell  volume  were  being 
incorporated  into  the  region.  Next,  a set  of  tracers  which 
coincide  with  the  interface  tracers  of  the  failed  material 
that  are  interior  to  the  cell  and  the  intercepts  of  the 
material  interface  with  the  cell  boundary  are  added  to  the 
tracer  string  and  the  exterior  cell  corner  points  are*deleted. 
Again,  the  interior  points,  the  intercepts  and ^ included  corner 
points  are  defined  by  CALFRC  and  VOLFND.  , r 


X Material  Trac«‘rs 
• Damaqc  Tracers 


1: 


STEP  I 

1.  Interior  damage  tracers  deleted. 

2.  Intercepts  of  damaqc  region 
boundary  and  cell  corner 
points  added  to  encompass 
entire  cell. 

STEP  II  , 1 

1.  Intercepts  of  material  inter- 

•’  face  and  duplicates  of  interior 
material  tracers  added. 

2.  Exterior  coll  corner  points 
deleted. 


I I 


I 

i 

I 
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The  reduction  of  a damage  region,  on  the  other  hand, 
is  accomplished  by  a single  call  to  CALFRC  and  VOLFND.  The 
damage  tracers  interior  to  the  failed  cell  are  deleted  and 
the  intercepts  and  included  corner  points,  as  defined  by 
CALFRC  and  VOLFND,  are  added.  \ 


The  steps  outlined  above  are  presented  in  greater 
detail  in  Figure  A. 2 which  is  an  illustrated  flow  diagram 
of  subroutine  MIXDMR.  , - 


NWDFLG  (called  from  VOLFND  and  MIXDMR) 


I 


-j 

When  the  boundary  of  a damage  region  first  intersects 
a cell,  subroutine  NWDFLG  is  called  to  assign  to  that  cell  a 
unique  value  of  DFLAG  and  to  initialize  the  locations  of  the 
IDDR  and  DVOL  arrays  that  will  be  associated  with  that  cell. 


The  DFLAG  array  functions  like  the  MFLAG  arrAy  by,,«^i|^- 
identifying  cells  that  are  exterior  to  all  damage  regions,  ' 
(DFLAG (k) »0) , cells  that  are  interior  to  N^^  damage  region 
(DFLAG(k) aN) , and  cells  that  contain  the  boundary  of  one  or 
more  damage  regions  (DFLAG (k) >1000) . In  the  latter  instance, 
the  value  of  (DFLAG (k) -1000)  is  used  to  define  the  second 
index  of  the  IDDR  and  DVOL  arrays  which  specify  the  identity 
and  volume  of  the  damage  regions  that  intersect  the  k^^’cell. 


NWDFLG  is  called  from  VOLFND  when  a daimage  region 
boundary  moves  -nto  a cell  that  was  previously  interior  to 
a single  region  or  was  outside  all  damage  regions.  NWDFLG  is 
also  called  from  MIXDMR  when  failure  first  occurs  in ^an  inter 
face  cell  and  a damage  region  is  created  whose  boundary  coin- 
cides with  a material  interface  that  intersects  the  cell. 


NWTCR  (called  from  PURDMR) 

When  a damage  region  is  created  from  a pure  cell  that 
was  partially  contained  in  one  or  more  damage  regions  that 
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are  being  reduced,  the  intercepts  of  those  reduced  regions 
with  the  cell's  boundary  are  saved  in  the  INTRCP  array  and 
added  to  the  tracer  string  of  the  new  region  by  subroutine 
NWTCR  > 


PURDMR  (called  from  SCALC) 

When  the  material  in  a pure  cell  that  is  hot  interior 
to  an  existing  damage  region  satisfies  the  NAG-FRAG  failure 
criterion,  subroutine  PUROMR  is  called  either  to  create  a 
new  damage  region  or  to  enlarge  an  existing  region.  In 
either  case  one  or  more  existing  damage  regions  may  be  reduced 

If  the  cell  is  not  already  intersected  by  a -lamage 
region,  PURDMR  simply  places  tracers  for  the  new  region  at 
the  four  corners  of  the  cell  and  calls  DAMAGE  to  define  the 
stresses  and  damage  parameters  of  the  newly  failed  material. 

If  less  than  half  of  the  cell  is  already  contained  in 
one  or  more  damage  regions,  these  regions  are  reduced  and 
the  tracer  string  which  defines  the  new  damage  region  con- 
tains the  intercepts  of  the  reduced  regions  as  well  as  the 
four  corner  points  of  the  cell  (see  NWTCR) . 

If  more  than  half  of  the  cell  already  belongs  to  one 
or  more  existing  damage  regions,  then  the  region  which  en- 
compasses the  largest  part  of  the  cell  is  enlarged  to  include 
all  of  the  cell  and  the  other  regions,  if  there  are  any,  are 
reduced. 

Region  N is  reduced  by  CALFRC,  VOLFND  and  INTSAV  when 
NCALL  = -N  and  N * NLSAVE  (HLSAVE  is  the  number  of  the  region 
being  enlarged) . CALFRC  and  VOLFND  compute  the  intercepts 
of  the  N^“  region  with  the  cell  boundary  and  determine  which, 
if  any,  corner  points  lie  between  the  entering  and  existing 
sides.  INTSAV  adds  these  intercepts  and  corner  points  to  the 
string  of  tracers  of  the  reduced  region  while  it  deletes  the 
damage  tracers  that  were  interior  to  the  cell.  PURDMR  rede- 
fines the  mass  of  the  reduced  region. 
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To  enlarge  a region,  PUKDMR  sets  NCALL  > -NLSAVE  and 
calls  CALFRC.  The  intercepts  of  the  NLSAVE  damage  region 
with  cell  (ISV.JSV)  are  computed  by  CALFRC;  VOLFND  identified 
the  included  corner  points.  INTSAV  redefines  the  tracer  string 
of  the  NLSAVE  region  to  include  these  intercepts  and  corner 
points  and  to  exclude  any  damage  tracers  that  were  interior 
to  cell  (ISV,JSV). 

The  Illustrated  flow  diagram  of  subroutine  PUROMR  in 
Figure  A. 3 provides  the  detailed  steps  in  the  procedure 
described  above. 

SCALC  (called  from  HELP) 

Subroutine  SCALC  first  calls  DAllAGE  to  update  the 
stresses  of  the  Lagrangian  damage  region.  Thereafter  SCALC 
is  concerned  with  the  computation  of  the  Eulerian  cell 
pressures  and  devlatoric  stresses.  Once  the  cell  stresses 
are  updated  by  CDT  and  DVSTRS,  the  NAG-FRAG  failure  criterion 
is  applied  to  those  cells  containing  undamaged  material. 

When  the  stresses  of  one  of  these  cells  satisfies  the  failure 
criterion,  SCALC  calls  either  PURDMR  or  MIXDMR,  depending  on 
whether  the  newly  failed  material  is  in  a pure  or  interface 
("mixed”)  cell.  PUROMR  and  MIXDMR  in  turn  either  create 
new  damage  regions  or  enlarge  existing  ones  to  encompass  the 
newly  damaged  material. 

After  the  new  damage  regions  have  been  created  and, 
in  some  cases,  the  existing  ones  modified,  the  stresses  of 
the  damage  regions  are  mapped  baclc  into  the  Eulerian  grid  as 
described  in  Chapter  3. 

A flow  diagram  of  SCALC  appears  as  Figure  A. 4. 

TPHASE  (called  from  HELP) 

As  part  of  the  changes  resulting  from  transporting 
material  in  two  passes:  first  radial  then  axial,  or  vice 
versa,  subroutine  TPHASE  was  created  to  call  the  transport 
routines,  INFACE  and  TRNSPT  for  each  pass  and  to  set  the 
flag  which  indicates  the  direction  of  the  pass.  TPHASE 
alternates  the  order  of  the  passes  each  cycle:  on  odd  numbered 
cycles  the  first  pass  is  in  the  axial  direction;  whereas,  on 
even  numbered  cycles  the  first  pass  is  in  the  radial  direction. 

(NOTE:  the  TPHASE  subroutine  in  the  documented  code 

was  renamed  TRNSPT  in  the  NAG-FRAG  version.) 

TRNVOL  (called  from  TRNSPT) 

The  volume  of  material  to  be  transported  in  one  direc- 
tion (radial  or  axial)  across  the  boundary  between  two  pure 


cells  is  computed  by  subroutine  TRMVOL.  As  described  in 
Section  3.2.1  the  volume  is  based  on  the  velocity  at  three 
positions,  two  on  the  cell  edges  parallel  to  the  direction 
of  the  transport,  and  one  at  the  cell  center.  (See  Section 
3.2.1.) 

TRVOLI  (called  from  DMCALC) 

Subroutine  TRVOLI  performs  the  same  function  a sub- 
routine TRNVOL,  except  one  or  both  of  the  two  colls  is  an 
interface  cell.  (See  Section  3.2.1.) 


A.  3 MODIFICATION  OF  EXISTING  HELP  SUBROUTINES 

Majo;;  modifications  were  made  to  several  of  the  existing 
HELP  subroutines  in  order  to  make  available  information 
required  by  the  NAG-FRAG  failure  model  and  also  to  improve 
the  synchronization  of  mass  and  tracer  motion,  the  latter 
being  necessitated  by  the  mapping  of  the  Lagrangian  damage 
region  stresses  into  the  Eulerian  grid.  The  modified  sub- 
routines are  discussed  below  in  alphabetical  order. 

CALFRC  (called  from  INFACE,  MIXDMR,  PURDMR,  SCALC,  SETUPA) 

In  the  documented  version  of  HELP  subroutine  CALFRC 
is  called  only  from  SETUPA  to  generate  the  interface  cells. 
(CALFRC  together  with  VOLFND  computes  the  volume  of  a given 
material  in  an  interface  cell  as  well  as  the  fractional  area 
of  the  right  and  top  faces  of  the  cell  belonging  to  that 
’.aterial.)  ’ 

In  the  present  version  CALFRC  and  VOLFND  are  used 
extensively  during  the  calculation  as  well  as  in  the  prob- 
lem generation.  CALFRC  is  called  at  least  once  each  cycle 
from  INFACE  to  compute  the  partial  volumes  'and  fractional 
cell  face  areas  of  materials  in  interface  (multimaterial  and 
free  surface)  cells.  The  partial  volumes  are  based  on  the 
position  of  the  interface  at  the  end  of  the  time  step.  Thus 
these  volumes,  together  with  the  mass  of  each  material  at 
the  end  of  the  time  step  (post-transport  phase) , determine 
the  new  densities  of  the  materials  in  interface  cells. 
(Previously,  the  new  densities  of  materials  in  a multimaterial 
cell  were  determined  in  an  iteration  which  equilibrated  the 
materials'  pressures;  the  densities  of  materials  in  a free 
surface  cell  were  not  changed  from  the  original  donor  cell 
value  unless  the  volume  so  defined  exceeded  the  cell  volume.) 

Once  at  least  one  Lagrangian  damage  region  has  been 
created,  CALFRC  is  called  from  SCALC,  MIXDMR  and/or  PURDMR 
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to  compute  the  volume  of  damaged  material  in  cells  containing 
the  boundary  of  one  or  more  damage  regions  and  to  compute 
appropriate  intercepts  and  corner  points  for  creating, 
reducing  and  enlarging  the  boundaries  of  existing  damage 
regions.  The  desired  function  to  be  performed  by  CALFRC  is 
indicated  by  the  value  of  the  blank  cononon  variable,  NCALL, 
as  discussed  below.  ) 

NCALL  « 0 Compute  the  partial  volumes  and 

fractional  areas  of  materials  in 
multimaterial  cells  (INPACE,  SETUPA) . 


NCALL  >=  NLOR  Compute  the  partial  volumes  of 

damage  regions  in  cells  containing 
damage  region  boundaries.  (NLDR 
is  the  number  of  deunage  regions 
that  have  been  created.)  (SCALC) . 

NCALL  = -N  (NiMXDMR)  Compute  the  intercepts  of  the  N^^ 

damage  region  with  cell  ISV,JSV; 
the  N^"  damage  region  is  being 
either  reduced  or  enlarged.  (MXDMR 
is  the  maximum  number  of  damage 
regions  that  can  be  generated.) 
(MIXDMR,  PUROMR) . 


NCALL  * -(MXDMR+1)  Compute  the  intercepts  of  the 

NMSAVE  material  package  boundary 
t . with  cell  ISV,JSV.  Store  the 

intercepts  and  included  corner 
points  in  the  XTS,YTS  arrays. 

(NMSAVE  is  defined  in  SCALC) (MIXDMR) 


When  processing  the  damage  region  boundaries  CALFRC 
examines  the  damage  region  tracer  arrays,  XP,  YP,  instead  of 
the  material  interface  tracer  arrays,  TX,  TY.  As  indicated 
above,  the  value  of  NCALL  instructs  CALFRC  as  to  which  arrays 
to  reference.  i j 

Material  interface  tracer 
arrays  are  referenced. 

i. 

Damage  region  tracer  arrays 
are  referenced. 

When  CALFRC  is  called  to  compute  the  volume  of  damaged 
material  in  the  Eulerian  cells  intersected  by  damage  region 
boundaries  (NCALL  « NLDR) , the  DFLAG  array  is  redefined  to 
reflect  the  new  position  of  the  damage  region  boundaries. 

For  example,  the  DFLAG  of  those  cells  no  longer  intersected 


NCALL  » 0 

NCALL  = -(MXDMR+1) 


TX,TY 


NCALL  » tN  (0<N<MXDMr)]xP,YP 
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by  a damage  region  boundary  is  set  to  zero  or<to  an  integer, 

N, (<NLDR)  depending  on  whether  the  cell  is  currently  outside 
all  damage  regions,  or  inside  the  region. 

CALFRC  defines  NTST  and  NTND  which  are  passed  to  VOLFND, 
then  INTSAV  through  an  argument  list  when  intercepts  are 
added  to  the  tracer  string  of  a damage  region.  NTST  and  NTND 
identify  the  damage  tracers  after  and  before  which  the  inter- 
cepts are  to  be  added  to  the  string.  - 


NTND 


Before  the  string  is  altered,  however,  CALFRC  insures  that 
the  beginning  tracer  of  the  string  lies  outside  of  cell(ISV, 
JS^  which  is  either  being  added  to  or  subtracted  from  the 
damage  region. 

CDT  (called  from  SCALC) 

In  both  versions  of  HELP  subroutine  CDT  computes  the 
pressures  of  the  Eulerian  cells,  using  the  Tillotson  equation 
of  state  (see  subroutine  EQST) , and  the  time  step  for  the 
current  computational  cycle  using  a Courant  stability  con- 
dition wherein  the  sound  speed  is  estimated  from  the  updated 
cell  pressures.  The  difference  in  subroutine  CDT  in  the  two 
versions  lies  in  its  treatment  of  multimaterial  and  free 
surface  cells. 

In  the  NAG-FRAG  version  a pressure  for  each  constituent 
of  an  interface  cell  is  computed  independently.  The  mass 
density  of  each  constituent  is  derived  from  its  current  mass, 
and  from  the  location  of  its  interface  in  the  cell.  The 
interface,  together  with  the  boundaries  of  the  cell,  specifies 
the  volume  of  the  material  Vj^,  and  the  density,  Pj^  is  simply 


M./V. . The  pressure  associated  with  the  center  of  an  inter- 
fi^ce^cell  is  a volume  weighted  average  of  the  pressures  of 
the  cell's  constituents. 

This  procedure  is  very  different  from  that  employed 
in  the  documented  version  of  HELP  in  which  the  position  of  the 
interface  has  no  influence  on  the  density  assigned  to  the 
materials  in  a mixed  cell.  Instead  the  material  densities 
as  well  as  the  cell  pressure  are  uniquely  determined  in  an 
iteration  which  adjusts  the  densities  of  the  cell's  constituents 
until  the  pressures  of  the  materials  are  equilibrated. 

The  two  procedures  give  rise  to  equivalent  stress 
fields  w^'^n  the  materials  of  contiguous  packages  are  similar, 
e.g.,  metals,  soils  and  polymers.  In  many  instances,  however, 
when  a metal-air  or  soil-air  interface  is  modeled,  the  revised 
procedure,  which  does  not  attempt  to  equilibrate  the  pressures 
of  the  constituents  of  mixed  cells,  is  preferable.  In  the 
present  scope  of  work  the  revised  procedure  was  chosen  because 
it  tied  the  calculation  of  volumetric  strain  to  the  position 
of  the  tracer  particles  in  the  Eulerlan  cells  as  well  as  in 
the  Lagrangian  damage  regions. 

CMPRSN  (called  from  SPHASE) 

Subroutine  CMPRSN  calculates  weighting  factors  that 
are  applied  to  cell-centered  stresses  when  the  stresses  at 
cell  boundaries  are  being  determined  (in  SPHASE) . The 
weighting  factors  prevent  cells  having  relatively  small 
masses  from  being  overaccelerated  by  more  massive  neighbor 
cells.  If  this  situation  occurs  primarily  at  the  free  sur- 
face, it  sufficies  for  the  weighting  factors  to  be  "com- 
pression," i.e.,  the  gross  density  of  the  cell  divided  by 
the  reference  density  of  the  material.  However,  when  a 
calculation  involves  large  density  discontinuities  between 
material  packages  (e.g.,  metal-air  interfaces)  the  "com- 
pression" weighting  of  cell  stresses  does  not  adequately  pre- 
vent overacceleration  of  the  less  massive  cells  which  contain 
normal  density  air. 

CMPRSN  has  been  modified  so  that  the  user  can  specify 
the  gross  density  of  a cell  instead  of  its  "compression" 
to  be  the  weighting  factor  in  the  calculation  of  cell  boundary 
stresses.  An  input  variable,  DENSWT,  has  been  added  to  the 
Z-block.  When  DENSWT  >=  1.0,  CMPRSN  returns  the  gross  density 
of  the  cell  Instead  of  its  "compression." 

DMADJ  (called  from  INFACE) 

The  mass  transport  terms  computed  by  DMCALC  for  inter- 
face cells  are  adjusted  by  subroutine  DMADJ  to  prevent  a cell 


from  being  "overemptied"  of  a material  and  to  exactly  evacuate 
a material  whose  interface  has  just  left  a cell. 

Because  the  mass  transport  is  executed  in  two  passes, 
the  adjustment  of  the  mass  transport  terms  on  the  second 
pass  must  only  affect  the  terms  associated  with  the  direction 
(axial  or  radial)  of  the  second  pass.  Tests  have  been  added 
to  DMCALC  to  perform  the  appropriate  adjustments  depending 
on  the  pass  number  (NFS)  and  the  direction  of  the  pass  (ID)  . 

FLGSET  (called  from  INFACE) 

Subroutine  FLGSET,  which  senses  when  a material  inter- 
face leaves  a cell  and  sets  the  appropriate  flags  for  that 
material  to  be  evacuated  by  DMADJ  and  in  some  cases  for  the 
cell  to  be  again  classified  as  "pure,"  has  been  rewritten  so 
as  to  use  the  VOLM  array  to  sense  the  presence  or  absence 
of  a material  interface  in  a given  interface  cell.  The  VOLM 
array,  which  stores  the  partial  cell  volume  occupied  by  each 
material  in  an  interface  cell,  was  not  computed  in  the  docu- 
mented version  of  the  code,  and  FLGSET  relied  on  values  in 
the  FRACRT  and  FRACTP  arrays  to  determine  whether  or  not  an 
interface  still  cut  across  a given  coll.  By  using  the  VOLM 
array,  FLGSET  has  become  a simpler  and  shorter  subroutine. 

FRACS 


Subroutine  FRACS  has  been  removed.  (See  CALFRC,  INFACE) . 
HELP  (main  driver  of  the  HELP  program) 

The  main  driver  routine  of  the  HELP  code  has  been 
modified  as  indicated  in  Figure  A.l.  Briefly,  the  changes 
reflect:  (1)  HPHASE  has  been  combined  with  SPHASE;  (2)  a new 

subroutine,  SCALC,  controls  the  computation  of  the  pressures, 
deviator  stresses  of  the  Eulerian  cells  as  well  as  of  the  damage 
regions;  and  (3)  the  pressures  are  no  longer  recomputed 'when 
the  calculation  is  restarted. 

HPHASE 

Subroutine  HPHASE  has  been  removed.  The  effects  of 
cell  pressures  on  cell  energy  and  momentum  are  no  longer 
accounted  for  in  a separate  computational  phase  of  the  cycle. 

The  pressures  are  added  to  the  normal  deviator  stresses  in 
SPHASE  and  the  effects  of  the  total  stress  components  are 
computed  there.  This  change  arose  from  the  need  to  con- 
sider and  modify  both  the  pressures  and  the  stress  Jeviators 
of  damaged  material  before  either  was  used  in  the  momentum 
and  energy  equations. 


INFACE  (called  from  TPHASE) 


Subroutine  INFACE  is  called  twice  each  cycle  from 
TPHASE,  once  for  the  radial  pacs  and  once  for  the  axial  pass. 
Because  the  transport  phase  has  been  split  into  two  directional 
passes,  it  is  less  necessary  to  subcycle  INFACE.  Furthermore, 
it  is  no  longer  possible  to  subcycle  INFACE  without  subcycling 
the  entire  transport  phase  of  the  cycle.  Thus  the  subcycling 
loop  has  been  removed  from  INFACE  although  the  subcycling 
counters  and  flags  have  been  retained. 

The  fractional  T'^olumes  and  cell  face  areas  associated 
with  the  materials  in  interface  cells  are  computed  by  CALFRC 
and  VOLFND  instead  of  by  FRACS  (subroutine  FRACS  has  been 
removed  from  the  program) . INFACE  calls  CALFRC  after  the 
transport  terms  have  been  computed  (DMCALC)  and  aftet  the 
tracers  have  been  moved  (MOVTCR) . Thus  the  fractional  areas 
used  to  compute  the  transport  terms  are  computed  on  the 
previous  call  to  INFACE.  On  cycle  zero  the  fractional  areas 
and  volumes  are  computed  as  the  mixed  cells  are  generated, 
i.e.,  when  CALFRC  is  called  from  subroutine  SETUPA. 

INPUT  (called  from  HELP) 

The  print  statements  and  tape  read  statements  in  sub- 
routine INPUT  have  been  extended  to  include  the  damage  param- 
eters of  the  materials  being  modeled  as  well  as  the  character- 
istics of  the  Lagrangian  damage  regions  that  have  been  created. 

The  mechanisms  for  sliplines  and  automatic  grid  rezoning 
have  temporarily  been  deactivated  by  statements  in  INPUT  which 
override  the  input  values  of  the  Z-block  variables,  NOSLIP, 

REZ,  NUMREZ,  lEXTX,  and  JEXTY. 

MOVTCR  (called  from  INFACE) 

To  better  synchronize  the  tracer  particle  motion  with 
the  mass  transport,  the  density  weighting  of  cell  velocities 
to  arrive  at  tracer  particle  velocities  was  removed  from  sub- 
routine MOVTCR.  In  addition,  the  logic  of  MOVTCR  was  altered 
slightly  to  accommodate  the  two  directional  subphases  of  the 
transport  phase:  MOVTCR  is  called  twice  each  cycle  and  moves 
the  tracers  in  only  one  direction. 

SETUP  (called  from  INPUT) 

Subroutine  SETUP  reads  input  cards  that  define  the 
TSR  array  for  each  material  in  the  problem  to  which  the  NAG- 
FRAG  failure  model  is  to  be  applied.  The  parameters  which 
govern  the  nucleation  and  growth  of  craclcs  (or  voids)  are 
stored  in  the  TSR  array.  Only  those  materials  whose  code 


numbers  are  1 through  6 can  utilize  the  NAG-FRAG  routines 
because  the  first  index  of  the  TSR  array,  which  identifies 
the  material,  is  dimensioned  six.  ‘ 

. ■ j. 

SPHASE  (called  from  HELP)  ' ^ ^ ' 

The  modified  SPHASE  subroutine  is  essentially  the 
second  half  of  the  same  routine  in  the  documented  version  of 
HELP.  The  first  half  has  been  replaced  by  subroutine  DVSTRS. 

Two  additional  options  have  been  added  to  SPHASE. 

First,  the  user  can  specify,,  through  a Z-block  variable, 

DENSWT,  whether  the  cell  boundary  stresses  will  be  a density- 
weighted  or  "compression"-weighted  average  of  the  cell- 
centered  stresses.  Second,  'the  user  can  activate  the  use  of 
a linear  and/or  quadratic  artificial  viscosity  when  it  is 
necessary  to  stabilize  the  stress  field.  The  artificial 
viscosity  terms  are  added  to  the  cell  boundary  stresses. 

In  the  present  version,  SPHASE  computes  the  effects 
of  the  total  stresses  (pressure  plus  stress  deviators)  on 
cell  momentum  and  energy;  whereas,  in  the  documented  version 
SPHASE  accounted  only  for  the  effects  of  the  deviator  stresses 


STRNG  (called  from  DVSTRS) 

Subroutine  STRNG  is  called  from  DVSTRS  to  compute  the 
shear  yield  strength  of  a given  pure  cell  or  of  a given 
material  in  a mixed  cell.  Formerly,  STRNG  computed  a volume 
weighted  average  yield  strength  for  mixed  cells,  but  this  is 
no  longer  necessary  in  the  current^ vei^kioh,  since  the  deviator 
stresses  of  each  material  in  a mixed'^cell  are  computed 
independently.  Similarly,  STRNG  no  longer  sets  to  zero  the 
yield  strength  of  materials  in  free  surface  cells.  The  void 
is  accounted  for  in  defining  the  cell-centered  stress  by 
volume  weighting  the  stress  deviators  of  the  constituents, 
including  the  void.  * 

\ 

TPHASE  (called  from  HELP) 

The  subroutine  called  TPHASE  in  the  documented  version 
of  HELP  has  been  renamed  TRN5PT.  The  new  TPHASE  routine  is 
simply  a driver  routine  which  controls  the  two  directional 
passes  of  the  transport  phase. 

VOLFND  (called  from  CALFRC) 

Normally  VOLFND  computes  the  fractional  cell  face 
areas  and  fractional  volume  of  a given  material  in  a given 
interface  cell.  VOLFND  takes  the  tracer  particles  that  are 


interior  ito  the  cell  and  the  intercepts  of  the  interface  h 
with  the  cell's  boundaries  (these  are  stored  in  the  XT,  YT  h 
arrays  by  CALFRC)  and  determines  which  corners  of  the  cell 
are  interior  to  the  specified  material  package,  thus  creating 
a set  of  points  which  circumscribe  the  fpart  of  the  material i ^ 
package  that  lies  inside  the  cell  in  question.  ” — " 

.)fl  I f 


Material  Package 


Define  Intercepts 
Interior  Points 


Add  jlncludocl 
Corners.  Com- 
I pute  fractional 
I cell  face  areas 
and  volumes. 


L 


».  t-  ^ 

In  , the  revised  version  of  HELP  these  mechanisms  are 
utilized  by  MIXDMR  and  PU^MR  to  modify  existing  damage  , ^ 

regions  to  include  or  exclude  portions  of  cells  that  have, 
failed.  They  are  also  used  to  compute  the  volume  of  damaged 
matarial  in  those  cells  containing  the  boundary,  of  one  or 
more  damage  regions. 

When  VOLFND  (which  is  always  called  from  CALFRC)  is 
being  used  to  enlarge  or  reduce  an  existing  damage  region,  it 
bypasses  the  calculation  of  fractional  areas  and  volumes. 

After  identifying  the  appropriate  cell  corner  points  and 
adding  them  to  the  XT,  YT  arrays,  VOLFND  calls  subroutine 
INTSAV.  Which  corner  points  are  chosen  in  these  cases  depends 
upon  whether  a region  is  being  reduced  or  enlarged,  ^l^hen  the^^ 
damage  region  is  being  reduced  the  array  of  interior^ points 
and  intercepts  provided  by  CALFRC  is  unaltered  and  the  same 
corner  points  are  chosen  by  VOLFND  as  would  be  chosen  if  the 
points  belonged  to  a materia  interface  rather  than  to  a 
damage  region  boundary.  On  the  other  hand,  when  a region 
is  being  enlarged  VOLFND  reverses  the  order  of  the  points 
provided  by  CALFRC  and  thereby  chooses  the  other  set  of  , 

corner  points. 
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The  intercepts  defined  by  CALFRC  and  the  corner  points 
defined  by  VOLFND  are  added  to  the ^damage  rdgibh  tracer  strings 
by  subroutine  INTSAV.  ^ \ v i \ - s f)n.  OC) 

. • iiirr  - u ' ■ • > • !-f!  . --  o 

' r ""  • 'I  . 70 


A.  4 DICTIONARY  OF  NEW  VARIABLES 

A number  of  new  blank  common  variables  and  new  arrays 
have  been  added  to  HELP  in  the  process  of  incorporating  the 
NAG-FRAG  subroutines  for  brittle  and  ductile  failure.  These 
new  variables  and  arrays  are  defined  below.  The  conventions 
used  in  describing  the  storage  location  of  the  variables  is 


The  variable  is  local  to  subroutine  NAME. 

The  variable  is  in  common  block  NAME. 

The  variable  is  in  Blank  Common  or 
equivalenced  to  a variable  in  Blank  Common. 

The  variable  is  equivalenced  to  a member  of 
the  Z-array,  the  first  array  in  Blank  Common. 
Most  of  these  variables  axe  used  in  generating 
and  restarting  problems. 

The  variable  is  used  as  a calling  argument. 

In  describing  the  dimensions  of  the  variables  the  following 
conventions  are  used: 


as  follows: 

(NAME) 

NAME 

B.C. 

= Z(N) 
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yjr 


The  variable  ia  not^dtroensioaed.’^  , 

< D •’t  * J i *f.  • ^ ^ 'jI 

The  array  is  always  dimensioned  30.  i 

The  dimension  of  the  array  mast  be  greater 
than  or  equal  to  the  value  of  NVAR#  where 
NVAR  is  determined  by  the  user  when  the 
calculation  ie' generated. 


Variable 


Dimension 


■O' 


Nc.tne 

Location 

of  Array 

Units 

DefiriltioVi 

CL 

■? 

^ PDMGE 

(5,MXDMR) 

3 

cm 

Cube  of  crack^^ifadius 

J 

parameter  for  each  of 

•'  t ' 

>T  r-per  r 

5 bin  angles. 

• V 

"1  : va  J 

4 

CN 

PDMGE 

(5,MXDMR) 

3 

cm 

Number  of  cracks  per 
cm^  for  each  of  5 bin 

'i  I > ‘j 

£‘-  ) 

I i 

** 

angles.  .so  ,'A'l,ov 

CVISR 

Z(126) 

Flags  for  generating 

CVIST  , 

Z(127) 

i-  f ■■  ■ : 

— 

reflective-grid  boundaries 
at  right  and  top,  respec- 

r ■• 

1 

tively. 

DENSWT 

Z(U) 

Flag  for  using  density 

1 : ■ ^ 

weighting  rather  than 

"compression"  weighting 

0*10  : 

when  defining  averaging  cell 

d a .J 

boundary  stresses. 

DFLAG 

.1 ' 

B.C, 

(KMAX) 

Flag  associated  with  each 

la  ■ 

i^A 

Eulerian  cell  which  indicates 

whether  the  cell  is  undamaged 

i i 

V i 

■ J : : ^ 

(DFLAG(K) =0) , is  interior 

T 

' t 

' ' ' 

to  a single  damage  region 
(DFLAG(K) <NLDR) , or  con- 

i t 

- 1 ■ 

: r^i 

tains  the  boundary  of  one 
or  more  damaged  regions 
(NLDR>DFLAG>100p;)^.  j i 

DRENG  , 

PDMGE  ; 

(MXDMR) 

ergs/g 

The  specific  internal  energy 

t 

'it 

of  the  material  in  each 

i _ i ' 

damage  region  is  stored 
in  the  DRENG  array. 

' a - 

) 

n ^v  ! TH30  '-i 

DRMASS 

, PDMGE 

(MXDMR) 

g 

The  mass  of  each  damage 

* ■ - 1 T . 

region  is  stored  in  the 

- V f'- 

DRMASJ  array. 

DS 

B.c! 

(IJDX1,2) 

cm 

Equivalenced  to  DX  and  DY. 

• ' k r 

IJDX1= 

MAXdNAX, 

Used  in  the  double  pass 
transport  phase*  i| 

4 

JMAX) +1 

tCJ*  - 

DSTLDR 

PDMGE  T xp 

' (3, MXDMR) 

dynes/ 

The  stress  deviators  of  the 

1 ' 

1 

cra2 

damage  regions  are  stored 

■1 

in  the  DSTLDR  array. 

Vi^rlable 

Name 


DVOL 


FARRAY 


fi 


FMODEL 


FRACBT 


IBGN 

TEND 

■ I 


Diitiension 

Iiocation  of  Array  Units  t.  A Definition  ©m  <4 


PDMGE  ,,  n (NLRMX, 
.V  NPDCLS) 

. : ' ) . n ■ 


^ 

PDMGE  (NMAT) 


■1  i 

--  ! V,  ‘ " 

p-  : . =i--  -..iw 


■ t-  V .f: 

‘ IM  '!  J ^ 

--  m :;a. - , 

, '1  , : : . I ' , ■ - t 

'.A.i-’u 

. - .2a 

T»>!-  T 

za28)  ^ — 


3 

cm  The  yolume  of  material  in  j'j 

an  Euieri >n  cell  belonging 
to  a given  damage  region  is 
stored  in  the  DVOr  array. 

(«i  (See  IDDR.)  aD.’  J ‘AO 


The  failure  model  to  be 
applied  to  each  material 
package  in  the  problem  is  IVD 
specified  by  the  user-assigned 
value  of  the  FARRAY  for  each' 
package;  e.g.,  for  package  N: 


FARRAY (N)  = 0. 

4 

FARRAY (N)  * 1. 

FARRAY (N)  = 2. 


Use  simple 
volumetric  * 
failure  cri- 
terion and 
zero  out 
stresses  in 
failed  cells. 

i 

Use  DFRACT 
(ductile 
material) . 

Use  BFRACT 
(brittle 
material) . 


An  input  flag;  when  non- 
zero it  indicates  that  one 
of  th6'  NAG-FRAG  ' failure  “ 
models  is  to  be  applied  to 
at  least  one  material  package. 


MXCELL  (NVOID, 

IMAX)'*' 

^ - 

J 

PDMGE  (MXDMR) 

PDMGE  (MXDMR) 


cm  The  fractional  area  sub- 

tended by  each  material 
package  interface  on  the 
bottom  boundary  of  each 
cell  in  the  bottom  row  of 
the  grid. 

!>■,  X 

The  Iseginning  and  ending 

index  of  the  tracer  co- 
ordinates (XP,YP)  associated 
,>iwith  a given  damage  region.^/’ 
'If,  are  stored  in  IBGN  and  lEND, 

respectively. 


Variable  Dimension 

Name  Location  of  Array 


Units 


Definition* 


IDDR 


PDMGE 


(NLRMX, 

NPDCLS) 


IGRP  t B,C. 


INTRCP  PDMGE 


(4,10) 


The  variable  'that  indicates 
which  direction  the  nass ,is 
being  transported; 

ID  * 1 radial 

ID  = 2 axial 

The  units  on  the  DO-loops 
in  TRNSPT  which  reflect  the 
direction  the  mass  is  being 
transported . 

The  complement  of  ID.  When 
ID=1,  then  ID3=2  and  vice 
versa.  ">  '*  !(U 

The  IDDR  array  identified 
which  damage  regions  inter- 
sect those  cells  that  con- 
tain one  or  more  damage 
region  boundaries,  e.g., 
if  MD  = DFLAG{K)-1000, 
then  IDD  = IDDR{N,MD)  N=l,2... 
and‘VL  = DVOLiN,MD)  NLRMX, 
where  VL  is  the  voi^ame  of 
cell  K that  is  inside  the 
damage  region  numbered  IDD. 

When  a given  material  inter- 
face crosses  a cell  more  than 
once,  two  or  more  non- 
contiguous regions  inside 
the  cell  belong  to  that 
material.  The  number  of^  ‘ 
such  regions  is  determined 
Ly  INTSAV  when  a given 
material  in  an  interface 
cell  fails.  When  1GRP>1 
more  than  one  damage  region 
must  be  created  to  encompass 
the  failed  material. 

The  intercepts  of  a damage 
region  with  the  boundaries 
of  a cell  are  stored  by 
INTSAV  when  the  region  is 
being  reduced. 


; ; : Kl  /id  = - fiV 

Variable  ; Dimension  to  ur  - btri  ' 

Name  Location  of  Arroy  Units  - ' — Definition — “ 


LBG  Z(106) 

LND  ZU07)^  j— - 

M • - (Y;  ‘’!  .■  . . -,-r 

. ' i:  - ; '1  liw  . " X 

' i(  . jS  t!  _ ' I ■ : ■ ' 

* • 1 


: X i 

LDRBGN  B.C.  

LDREND  B.C.  

_ ^ I J ^ . 

: p xf 

ir  ' - r I ' I ^ 

; r : ' ; C 1 0 i 

, . • - It;  J y.l 

, tv- ; ~(  y ij 

I’l,  .)  1 

LSLDR  (CIM  B.C.  

\y 

. t -j  :■  I V 

• J - - . V ' . 


NBSAVE  B.C. 

NESAVE  B.C. 


NCALL  ~L  B.C.r  r n^x--- • 

■J-  1 '■  ’< 

lix  V ■.-■Vit 

n i 


The  indices  of  the  tracers 
of  a given  damage  region 
that  just  precede  (LBG) 
and  just  follow  (LND)  those 
tracers  that  are  inside 
cell  (ISV.JSV)  . When  a'r^l 
damage  region  boundary  is 
being  enlarged  or  reduced, 
LBG  and  LND  identify  where, 
in  the  XP,  YP  arrays,  tracers 
should  be  added  or  deleted. 

Subroutine  DAMAGE  is  called 
either  to  update  all  of  the 
. dimage  regions,  in  whichH 
OJease  LDRBGN»1  and  LDREND* 
*NLDR,  or  to  establish  the 
stresses,  etc.,  of  a new 
damage  region,  in  which 
case  LDRBGN=LDREND. 

The  initialization  of  the 
arrays  used  by  BFRACT  occurs 
when  LSLDR*0,  i.e.,  when  a 
calculation  is  generated 
or  restarted.  LSLDR  is  set 
to  zero  in  INPUT  and  reset 
to  one  in  BFRACT. 

These  are  the  indices  of  the 
entering  and  exiting  inter- 
cepts of  a given  damage 
region  with  cell  (ISV,JSV) 
after  the  region  lioundary 
has  been  enlarged  or  reduced 
to  include  or  exclude  cell 
(ISV,J&V)  and  after  the 
appropriate  intercepts  and 
corner  points  have  been 
added  by  INTSAV. 

A flag  defined  by  SETUP, 
MIXDMR,  PURDMR  and  SCALC 
which  establishes  what 
function  CALFRC  is  to 
perform  (see  discussions 
of  CALFRC  in  A. 2) . 


100 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


V 


I 

j 


NEXIT  B.C. 

> fL  1 

' I vii  "ji  ^ ( r 

•■I  ^ ; .>1 


i - i -riiTi 

i;  ;■-.  _ ■ i 


NEXT  ^ PIMGE  5 

? n ■ : ; 1 


NLDR  Z<147) 


NLSAVE  B.C. 


NMPDR  (MXDMR) 


-rj 


■\  i.:  -iii  ''  < 

Next-is  the  index-of  the 
exiting  intercept  in  ,the  . 
XTS,YTS  arrays'of  a given' 
material  interface  with 
cell  (ISV,JSV).  (Equals 
the  number  of  points  (NTEMP) 
transferred  by  CALFRA  to 
VOLFND.  Used  in  INTSAV  to 
define  the  NEXT  array. 


.1  t 

The  index  of  the  exiting 
intercept  of  each  region 
of  the  NMSAVE  material  in 
cel?  (ISV,  JSV)  is  stored 
in  the  NEXT  array.  (See 
NEXIT.) 


The  number  of  existing 
damage  reyions.  Incre- 
mented ii.  MIXDMR  and  PUI^MP. 

' , (initially  zero)  . ' 

The  number  of  the  damage 
region  being  enlarged  to 
, anepmpass  all  of  cell  (ISV, 
JSV)  (PURDMR)  or  all  of  the 
NMSAVE  material  pac)cage  in 
cell  (ISV, JSV)  (MIXDMR) . 

The  material  )pacl(age  number 
of  each  damage  region  is 
stored  in  th6)  NMPDR  array. 
From  this  variable  the  code 
number  of  the  material  of  a 
damage  region  is  found,  e.g., 

N = NMPDR(I) 

'MN  = MAT(N)  Vfi  , V >! 

MN  is  the  code  number  of  the 
material  in  the  I^“  damage 
region. 

The  material  paclcage  number 
of  the  material  in  a mixed 
^ cell  that  has , satisfied  the 
NAG-FRAG  failure  criterion. 
Defined  in  SCALC;  used  in 
MIXDMR. 


-•  .'ySx 


f 


•I 


V r 


(^O  I i 


Variable  ' ’ '*  Dimension^ 

Name  Location  of  Array 


if  .1 

tU  yr;  1 1 A 

Units 


O for- 1 1 
i i£,H 

DefinitJon 


NPDCLg  i Z(X36V 

■•I  T i i. 


' iw  z>i  ■■  1^.-  f 

. {V  • .■ 

KC  i':- 1 V i o<{  : ui  ..  i 

■hmjad  ya  V - C'*’ 

; f V f**I  r!  ’ -s  , 

. .'tTfifi  ^y'iH  /rt  * " > f 


NTS 


I ft  ■ 


i f'( . 


PRS 


1 111 


■I  ■!■) 

^ ; 


PRSLDR  ’ 


■V 

r*.-;  V - » 


QLIN 

QQUAO 


RLVSAV 


ri 


ROT 


fivK'i 


RVOL 

1 , ,1 


vJ  1/  '.^OP  . 


k \ ’ p V : t;: 


e*’  t i 
' AnM"!  >f\  ! 

• L , • . . - 

. r”j;;  ^ 

f •:  : 


MXCELL  (NMAT, 


NMXCLS) 


PDMGE  " (MXDMR) 

1 


2(93) 


( 


2(94) 


PDMGE 


(MXDMR) 


PDMGE 


(MXDMR) 


; i 

i O ” . A ■ 


PDMGE  ‘1  (MXDMR) 

^ I - . 


The  maximum  number  of  cells 
containing  one  or  more 
damage  region  boundaries 
on  any  one  cycle.  NPDCLS 
should  correspond  to  the 
second  dimension  of  the 
DVOL  and  IDDR  arrays. 


The^ number  of  material 
tracers,  intercepts  and'^^^'^ 
corner  points  (defined  by 
VOLFND  and  stored  in  the 
XTS,  YTS  arrays  by  INTSAV) 
which  circumscribe  each  of 
at  most  five  regions  of  the 
NMSAVE  materi-al  package  in 
cell  (ISVi-U6v). 


HC 


dynes/ 


cm^ 


The  pressure  of  each 
material  in  an  interface 

cell  is  stored  in  the  PRS 

.iK 


array. 


dynes/ 


cm 


The  pressure  of  each  damage 
region  (as  defined  by  BPRACT 
or  DFRACT)  is  stored  in  the 
PRSLDR  array. 


The  coefficients  of  t^ie^ 
linear  and  quadratic  arti- 
fic^.al  viscosity  terms, 
respectively,  added  to  the 
cell  boundary  stresses  in 
SPHASE  (see  Section  4.2). 


The  relative  volume  of  each 
damage  region  before  the 
BFRACT  or  DFRACT  are  called. 


radians  The  cumulative  rotation  of 

each  damage  region  is ^stored 
in  the  ROT  array.  ’ 


The  relative  volume  of  each 
damage  region  (updated  in 
DAMAGE)  based  on  the  region's 
gross  density. 
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Variable 

Name 


RWLDR 


SAMMXY 

SM4PXy 


SAMMUV 

SAMUVT 


STRR^ 

STRZ 

STZZ 


TSR{N,1) 


Dimension 

Location  of  Array 


nC-CJi 


Units 


Definition 


PDMGE 


B.C. 


MXCELL 
by  X 

equiva- 

lence 


MXCELL 

by 

equiva* 

•^lence 


MXCELL 


(1  w 

PDMGE 


(MXDMR) 


j The  relative  void  or  crack 

' J L J . 


volume  in  each  damage  region 
(updated  in  BFRACT  and  DFRACT) 


Maximum 

of 

(IMAX, 

JMAX)tl 

• /I 


Equivalenced  to  X and  Y. 

Used' in  double  pass  trans- 
port phase. 

.vt>Mdv 

(.''.WiRST 


(NMAT« 

NMXCLS, 

2) 


The  mass  of  each  material 
to  be  transported  across 
the  left  Slid  right  boundaries, 
respectively,  of  a mixed 
cel]  on  the  radial  pass  of 
the  transport  phase,  and 
across  the  bottom  and  top  


boundaries,  respectively. 


, . w , 

on  the  axial  pass.  Equiva-  ^ 
lenced  to  SAMMY,  SGAMC  and 
to  SAMPy,  SAMMP,  respectively. 


(Nf4AT, 
NMXCLS, 
2)  . 


g-cm/ 

sec 


The  momentum  of  each  material 
to  be  transported  across  the 
boundaries  of  a mixed  cell. 
Equivalenced  to’SAlMMU,  SAMMVT 
and ^ to  SAMUT,  SAMVT, 
respectively. 


(NMAT, 

NMXCLS) 


dynes/ 

r»m^ 


The  stress  deviators  (S  , 
^rz'  ^zz^  associated  wiiK 
each  material  in  a mixed 


each  material  in  a mixed 
cell. 


(6,30) 
Storc'ge 
available 
for  at 
most  six 
materials 
in  a cal-  , 
culation.  ® 


cm^/ 


dynes- 


Void  (or  crack)  growth  con- 
stant = 4/3n;  n is  material 
viscosity. 


T r» . • '• 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


TSR(N,2) 

\ 


TSR(N,3) 

. > 

TSR{N,4) 

through 

TSR(N,e) 


i ■ 1 

■ i.hi  - 

f-  , 


TSR(N,9) 

through 

TSR(Na3) 

. VI  .3V. 


TVSAVE  ^ 


^ 1 


UNFRAG 


UV 


- . J I <1  ' i Ki 

PDMGE  All  30  

locations 

each 

PDMGE ' ^ material 

^ are  riot 

- oj  b currently 

j oirlu  (j 

PDMGE 


■i  o 


'] 

I. ; 


'J  V - > i ^ . 1 ^ 

dyne/  r^.Void  (or  crack)  growth 
cm^  threshold. 

cm  Void  (or  crack)  nucle- 

ation  radius  parameter. 


int 


no,/ 

cm-^-sec  . 

dynes/  *^no'°no' 
cm^  . . 

dnes/  ^1^°1^ 
cm^ 


Used  in  void 
nucleation  rate 
equation. 


N=N 


e(P-Pno)/Pi 


(Bfract) 


N»N  e^‘^'J>'l'"®no)/Pi 

(Sfract)  r^y[}AA’ 
where  is  the 
stress  normal  to 
the  plane  of  the 
cracks. 


PDMGE 

■ .j  .id 

. <?  Ir.  5/;; 

. I 

t ( '1A  J % Y 1 ■ i 


Fragmentation  parameters 
provided  by  SRI. 


. i>.  ;n  ''j.:c'T  criT 

•v;  ‘ 3 t oi;:  1 J 

i:  m 3 « .n  * 

-B.C.  Maximm  of  cm^ 

(IMAX.JMAX) 

Y +1 

' - VOU  Ifc  ! ._  ' \ ! 

■ . f'.  -i 

n ^ i ^ 1 if' 

. i 


PDMGE  (MXDMR) 


J,  VUfIlUr: 

■'■tO--  vd 

V u 

The  transport jVolume  com- 
puted between  pure  cells 
is  based  on  the  velocity 
at  three  positions.  The 
velocity  at  the  top  (ID*1) 
or''left  (ID=2)  is  saved  lin 
TVSAVE  for  use  when  the  Tc' 
next  row  (ID=1)  or  column 
(ID“2)  of  the  grid  is 
■processed,  d-  •]  (i,ii)M?.T 

The  fraction  of  unfragmented 
material  in  each  damage 
region  of  a brittle  material. 


B.C. 


n.f  - 

(KMAX,2)  cm/sec  Equlvalenced  to  U and  V. 

Used  in  the  double  pass 
transport  phase. 


f 
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y-j;  ■ i : : T ZV  ‘lV/V  i-i  : ::  A 

Variable.!  . . Dimenaion  ^ r ■ s i \ 

Name  L >catlpn  { of  Array  Units  Definition  ij-^,3 

: — -.-T-! : — ■ — . .'■  " 

VLFAIL  i 25  (98)  1 »*>o»r -i-.c  w._  The  Volume  of  unfailed  1 

» =fr(j  fo  r ■<:.!  j c material  must  Be  greater 

n f o io  r-n  than  VLFAIL  times  the  volume 


-x-j’O?  •A'j'T  1 i f Of  that  material  in  the  cell 


' , ■ 1 

:fa  before  a new  damaged  region 

is  created  or  an  existing 

. v;«'  1 

• fa->  :.  1 

region  is  enlarged. f (Default 

»IJ 

-- 

» A. 

value  * 0.0?.)  rt 
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NMXCLS) 

cm 

The  volume  of  each  material 
in  an  Interface  cell,  as  com- 

-5 

'^puted  from  the  positions  of 

- the  tracer  particles  in 
^CALFRC  and  VOLFND,  is  stored 

.( joo 

e.’.-a 
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l afl  ti.a  ■. ' 

centered  pressure  Of  a 
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■ , • i -- 
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an  interface  cell. 
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MODIFICATIONS  TO  THE  HELP  INPUT 


A. 5 


A minimum  number  of  changes  were* made  to  the  HELP  inputt^ 
Some  Z-block  V2lriables  were  added,'  some 'were  'deleted' and  some - 
were  relocated.  A set  of  cards  defining  the  NAG-FRAG  niicleation 
and  growth  pareuneters  for  each  material  being  modeled -by  either^  > 
BFRACT  or  DFRACT  was  added  to  the  input,  and  the  input  cards 
defining  the  shearAstrength  and  tensile  strength  of  each 
material  package  ’were  altered  to  define  another  array  speci- 
fying the'  failure  model  toibe  employed. 

The  Z-block  variables  that  were  added  are  listed  below. 
They  are  defined  in  the  pireceding  dictionary  and  in  the  input 
form  of  Appendix  B. 

' ^ /M)  . K .l'i 


DENSWT 

CVISR 

MXDMR 

CVIST 

QLIN  , 

-4  ')FMODEL 

QQUAD 

NPDCLS 

VLFAIL 

NLRMX 

VOLMIN 

a 

t. 

The  deleted  Z-block  variables  were  related  either  to  the 
removal  of  the  plugging  option  and  the  subcycling  of  iNFACE, 
or  to  the  modified tartificial  viscosity; 

r 

^ plugging  failure' removed  • i ' 


1 i 

I { JNFACE  not  subcycled 

)use  of  artificial  viscosity  now 
^ Specified  by  QLIN  and  QQUAD. 

The  variables  that  were  relocated  in  the  Z-block  so  as  to  be 
next  to  closely  related  parameters  were; 

CVIS  (moved  from  Z(27)  to  Z(125)  to  be 

near  CVISR  (Z(126))  and  CVIST  (Z(127)) 

CRATIO  (moved  from  Z(148)  to  «(61)  to  be 

near  DENSWT  (Z(60)) 

The  cards  that  define  the  shear  yield  and  tensile 
strength  of  each  material  package  now  define  an  additional 
array,  FARRAY,  which  specifies  which  failure  model  is  to  be 
applied  to  each  package.  (See  the  input  form  in  Appendix  B.) 


IPLGRT 

PLWMIN 

IPLGBT 

IPLGTP 

PLGOPT 

CYCMX 

■ ^ 

LVISC 
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The  TSR  variables  that  store  the,NAG-FRAG  nucleation 
and  growth  parameters  are  defined  by  the  fourth  set  of  input 
cards  following  the  Z-block  data.  They  are  defined  in  the 
dictionary  as  well  as  in  Appendix  B.  | 

This  version  of  HELP  can  be  run  without  utilizing 
either  BFRACT  or  DPRACT.  When  those  models  are  employed,  how” 
ever,  the  grid  cannot,  at  this  time,  be  automatically  rezoned. 
Thus  when  FMODEL  = 1.,  uWe  rezone  parameters  are  turned  off  in 
subroutine  INPUT.  * j 

■ i ’ I 

The  subroutine  which  approximated  the  slippage  between 
materials  separated  by  a slip  line  (UVCALC)  was  removed 
because  it  depended  upon  momentum  terms  resulting  only  from 
the  pressure  field  as  computed  by  HPHASE.  Because  the  effects 
of  the  pressure  and  stress  deviators  have  been  combined  in 
SPHASE,  the  simple  algorithm  employed  by  UVCALC  must  be  re- 
placed by  one  which  considers  the  components  of  the  total 
stress  normal  and  tangential  to  the  slip  line  in  the  cell. 
Until  this  modification  is  made  and  a subroutine  added  which 
replaces  UVCALC,  slip  lines  cannot  be  used.  Thus  subroutine 
INPUT  overrides  all  parameters  which  attempt  to  generate  and 
use  slip  lines.  r?’. 

■ lOV  , 

T-  , 

A.  6 STORAGE  REQUIREMENTS  OF  HELP  WITH  NAG-FRi;G  FAILURE  | 

MODELS 

When  segmented,  as  shown  in  Figure  A.  5,  and  executed 
on  the  UNI  VAC  1108,  the  version  of  the  HELP  code  with  the 
NAG-FRAG  failure  models  and  the  associated  systems  routines 
require  approximately  17,800  words  of  core  storage  for  the* 
instruction  bank.  This  represents  an  increase  of  about  6800 
words  over  the  1975  documented  version  of  HELP. 

In  the  figure,  the  subroutines  in  the  first  column, 
which  is  labeled  LEVEL  1,  are  loaded  into  core  and  remain  in 
core  for  the  duration  of  the  calculation.  The  blocks  of  sub- 
routines, or  segments,  under  LEVEL  2 and  LEVEL  3 are  then 
rolled  in  and  out  of  core  as  needed  by  the  various  stages  of 
the  computational  cycle.  The  INPUT-GENERATOR  segment  is 
executed  only  once  per  run,  either  in  generating  or  restarting 
a calculation.  In  most  of  the  cycles  only  three  LEVEL  2 seg- 
ments (the  SCALC  segment,  the  SPPASE  segment,  and  the  TPHASE 
segment) ! and  four  LEVEL  3 segments  are  rolled  in  and  out  of 
core.  (At  the  present  time  the  REZONE  routine  is  deactivated 
and  so  the  REZONE  segment  is  never  loaded.) 

The  amount  of  storage  required  for  the  data  bank,  in 
this  new  version  of  the  code,  is,  of  course,  highly  dependent 
on  the  dimensions  of  the  grid,  the  number  of  materials  and 
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LEVEL  2 LEVEL  3 

INPUT-GENERATOR  SEGMENT  i-  ‘ .il  dA 


CAROS 

jT  * 

SETUP A - 

INPUT 

FILGRD  ' 

SETUP 

TSETUP  , 

? -.r 

Figure  A. 5. 


'rll.  r/  V 

A schematic  representation  of  an  efficient  method  'o 
of  segmenting  HELP  which  minimizes  core  storage 
requirements. 
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the  extent  of  the  damage  incurred  in  the  course  of  the  calcu- 
lation. The' number  of  cells' that  at  one' time  or  another  will 
satisfy  the  NAG-FRAG  failure  condition  will  determine^  the 
number  of  damage  regions  that  will  be  created  (MXDMRK  and  the 
maximum  number  of  Eulerian  cells  that  will  contain' damage 
region  boundaries  on  any  one  cycle  (NPDCLS) . 

Twelve  arrays  which  store  information  on  the  damage 
regions  have  been  added  to  the  data  storage  requirement  of 
HELP.  If  the  user  estimates  that  at  most  100  damage  regions 
will  be  created  then  these  arrays  would  add  22  x 100  or  2,200 
words  to  the  data  bank.  In  addition  to  these  twelve  arrays, 
two  more  arrays,  DVOL  and  IDDR,  which  store  information  on 
cells  containing  damage  region  boundaries,  require  core  storage. 
Since  one  can  expect  that  no  more  than  nine  different  damage 
region  boundaries  will  cross  a given  cell  on  the  same  cycle, 
the  first  index  of  DVOL  and  IDDR  can  safely  be  nine.  There- 
fore, if  one  expects  at  most  100  cells  to  contain  damage 
region  boundaries,  then  another  2 x 9 x 100  or  1,800  words 
would  be  added  to  the  data  bank. 

In  addition  to  these  arrays  whose  dimensions  depend 
upon  the  volume  of  damaged  material,  there  are  five  arrays 
(INTRCP,  NEXT,  TSR,  XTS,  YTS)  whose  dimensions  are  fixed; 
these  arrays  require  an  additional  520  words  of  core  storage. 

The  DFLAG  array  was  added  to  blank  common  in  order  to 
identify  those  cells  containing  damaged  material,  those  con- 
taining both  damaged  and  undamaged  material,  and  those  con- 
taining only  undamaged  material.  The  DFLAG  array  is  dimensioned 
by  the  total  number  of  cells  in  the  grid.  In  addition,  four 
arrays  (PRS,  VOLM,  FRACTP,  FRACBT)  were  added  to  the  named 
common  block,  MXCELL.  Their  dimensions  depend  upon  the  number 
of  materials  and  the  maximum  number  of  interface  cells  that 
exist  on  a single  cycle  of  the  calculation.  The  PRS  and  VOLM 
arrays  store  the  pressure  and  volume  of  each  material  in  each 
interface  cell.  The  fractional  cell  area  arrays,  FRACRT  and 
FRACTP  were  previously  equivalenced  to  the  P array,  but  since 
they  must  now  be  saved  between  cycles,  they  too  have  been 
added  to  MXCELL. 

To  summarize  the  effects  on  the  storage  requirements 
of  the  revised  HELP  code,  brought  about  by  the  additional 
arrays  described  above,  assume  the  following  problem  parameters; 

1500  cells 

2 material  packages 

300  interface  cells 

100  damage  regions 

100  cells  containing  damage  region 

boundaries. 
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A calculation  with  these  pareuneters/ using 'the  version  of  " ^ ' 
HELP  with  the  NAG-FRAG  models,  would  require  (2200  + 1800 
+ 520  + 1500  +.11.  X 300)  * 9320  more^  words  in  the  data  bank 
than  woulci  the^  same  calculation  using  the  1975  documented  V _ 
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APPENDIX  B 


INPUT  FORM  FOR  HELP  WITH  NAG-FRAG  FAILURE  MODELS 


• The  Input  form  Included  here  reflects  the  modifications 
discussed  in  the  previous  section.*  The  default  values  of  most 
Z-block  variables  appear  in  parentheses  after  the  definition; 
if  no  value  is  specified,  the  default  value  is  zero  as  blank 
common  is ^initialized  before  the  input  cards  are  read. 
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CARDS  DEFINING  DX  AND  DY  ARRAYS 


CARDS  DEFINING  INITIAL  CONDITIONS  OF  EACH  MATERIAL  PACKAGE 
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DEFINING  NAG-PRAG  PARAMETERS  OF  EACH  MATERIAL  PACFAGE 


il 

% 

z 

% 

z 

» « 
z z 

% 

z 

z 

s « 

hi  *0 

ec 

cc 

» ee 

z 

z 

CARDS  DEFINING  TRACER  PARTICLES  OF  EACH  MATERIAL  PACKAGE  AND  THE  VOID 

Repeat  cards  1,  2 and  (if  required)  3 for  each  line  segment  or  arc  of  each  material 
and  void  package  boundary  (except  package  boundaries  that  are  also  grid  boundaries). 

for  each  package  must  be  input  in  an  order  that  will  put  the 
package  to  the  left  of  any  two  consecutive  particles.  Refer  to  Section  7.2.S. 
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CARDS  DEFINING  SLIPLINE  ENDPOINTS 
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ATl’N:  DRDMI-R 

DRDMI-RBL  ,, 
Redstone  Arsenal,  AL  35809 

I « 

1 Commander 

US  Army  Missile  Materiel 
Readiness  Command 
ATTN : DRSMI -AOM 
Redstone  Arsenal,  AL  35809 

1 Congnander 

US  Army  Tank  Automotive 
Research  8 Development  Cmd 
ATTN:  DRDTA-UL 
Warren,  MI  48090  j 


1 Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRSAV-E 
12th  and  Spruce  Streets 
Sr.  Louis,  MO  63166 

i Director 

US  Array  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1 Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  .NJ  07703 


2 Commander  !- 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-TSS  (2  cys) 
Dover,  NJ  07801 
. . 

4 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  Mr.  V.  Guadagno 

Mr.  R.  Davitt  , 
Dr.  J.  T.  Frasier 
G.  Demitrak 
Dover,  NJ  07801 

] Commander 

US  Array  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L,  Tech  Lib 
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5 Commander 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN:  DRXMR-T,  Mr.  J.  Bluhm 
Dr.  D.  Roylance 
Dr.  A.  F.  Wilde 
Dr.  J.  Mescal 1 
DRXMR-ATL 

Watertown,  MA  02172 

1  Director  > 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range 
NM  88002 

1 Deputy  Assistant  Secretary 
of  the  Army  {R8D) 
Department  of  the  Array 
Washington,  DC  20310 

1 HQDA  (DAMA-ARP) 

Washington,  DC  20310 

1 HQDA  (DAMA-MS) 

Washington,  DC  20310 

1 Commander 

US  Army  Research  Office 
ATTN:  Dr.  E.  Saibel 
P.  0.  Box  12211 
Research  Triangle  Park 
NC  27709 


1 Office  of  Naval  Research  I 
ATTN:  Code  ONR  439,  N. Per rone 
Department  of  the  Navy 

800  North  Quincy  Street 
Arlington,  VA  22217 

3 Commander  i 

Naval  Air  Systems  Command 
ATTN:  AIR -604 
Washington,  DC  20360 

3 Commander 

Naval  Ordnance  Systems  Command 
ATTN:  ORD-9132 
Washington,  DC  20360  i 

2 Commander  i , 

Naval  Air  Development 

Center,  Johnsville 
Warminster,  PA  18974 

( 

1 Commander 

Naval  Missile  Center 
Point  Mugu,  CA  93041 

1 Commander  8 Director 

David  W.  Taylor  Naval  Ship 
Research  8 Development  Ctr 
Bethesda,  MD  20084 

1 Command<»r 

Naval  I urface  Weapons  Center 
ATTN;  Code  TX,  Dr.  W.G.  Soper 
Dahlgren,  VA  22448 


1 


Commander 

US  Army  BMD  Advanced 
Technology  Center 
ATTN:  BMDATC-M,  Mr.  P.  Boyd 
P.'  0.  Box  1500 
Huntsville,  AL  35807 


f 


2 Cononander 

Naval  Surface  Weapons  Center 
Silver  Spring.  MD  20910 

3 Commander 

Naval  Weapons  Center 
ATTN:  Code  4057 

Code  4011,  Dr.E.Lundstrom 
Code  3813,  Mr.M.Backman 
China  Lake,  CA  93555 


I 


I 

« 4 


i 

o ! 
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4 Commander 

Naval  Research  Laboratory 
ATTN:  Mr.  W.  J.  Ferguson 
Dr.  C.  Sanday 
Dr.  H.  Pusey 
Dr.  F.  Rosenthal 
Washington,  DC  20375 

1 Superintendent 
Naval  Postgraduate  School 
ATTN:  Dir  of  Lib 
Monterey,  CA  93940 

2 ADTC/DLJW  (MAJ  D.  Matuska 

LTC  J.  Osborn)  ‘ 

Eglin  AFB,  FL  32542 

1 AFFDL  (FDT) 

Wright-Pattf rron  AFB,  Ovi  45433 

1 AFML  (Dr.  T.  Nicholas) 
Wright-Patterson  AFB,  OH  45433 

2 ASD  (XROT,  Gerald  Bennett; 

ENFTV,  Martin  Lentz) 
Wright-Patterson  AFB,  OH  45433 

3 Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  R.  H.  Toland,  L-424 
Dr.  M.  L.  Wilkins 
Dr.  R.  Werne 
Livermore,  CA  94550 

1 Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 

1 Director 

Jet  Propulsion  Laboratory 
ATTN:  Lib  (TD) 

4800  Oak  Grove  Drive 
Pasadena,  CA  91103 
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4 Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
‘Langley  Station  ‘ 

Hampton,  VA  23365 

1 Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
ATTN : Lib 
Houston,  TX  77058 

i 

1 Aeronautical  Research 
Associates  of  Princeton,  Inc. 

50  Washington  Road 
Princ  ‘ton,  NJ  08540 

2 Aerospace  Corporation 
ATTN:  Mr.  L.  Rubin 

Mr.  L,  G.  King 
2350  E.  El  Segundo  Blvd 
El  Segundo,  CA  90009 

1 Boeing  Aerospace  Company 
ATTN:  Mr.  R.  G.  Blaisdell 
(M.S.  40-25) 

Seattle,  WA  98124 

1 Effects  Technology,  Inc. 

5383  Hollister  Avenue 
P.  0.  Box  30400 
Santa  Barbara,  CA  93105 

1 Falcon  R^D 

ATFN:  Mr.  R.  Miller 
1225  S.  Huron  Street 
Denver,  C9  80223 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1  FMC  Corporation 

Ordnance  Engineering  Div. 

San  Jose,  CA  95114 

1 General  Electric  Company 
Armament  Systems  Dept. 
Burlington,  VT  05401 

1 President 

General  Research  Corporation 
ATTN : Lib 
McLean,  VA  22101 

1 Goodier  Aerospace  Corp. 

1210  Massillon  Road 
Akron,  OH  44315 

3 Honeywell,  Inc. 

Government  § Aerospace 
Products  Division 
ATTN:  Mr.  J.  Blackburn 
Dr.  G.  Johnson 
Mr.  R.  Simpson 
600  Second  Street,  NE 
Hopkins,  KW  55343 

1 H.  P.  White  Laboratory 
Bel  Air,  MD  21014 

1 Kaman  Sciences  Corporation 
ATTN;  Dr. ,P.  Snow 
1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80933 

I 

1 Lockheed  Corporation 
ATTN:  Dr.  C.  E.  Vivian 
Department  8114 
Sunnyvale,  CA,  S4087  , 

1 Materials  Research 
Laboratory,  Inc . 

1 Science  Road 
Glenwood,  IL  60427 
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1 McDonnell -Douglas 

Astronautics  Company 
ATTN:  Mail  Station  21-2 
Dr.  J.  Wall 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

1 New  Mexico  Institute  of 
Mining  and  Technology 
Terra  Group 
Socorro,  NM  87801 

1 Pacific  Technical  Corporation 
ATTN:  Dr.  F.  K.  Feldmann 
460  Ward  Drive 

Santa  Barbara,  CA  93105 

2 Physics  International  Company 
ATTN:  Dr.  D.  Orphal 

Dr.  E.  T.  Moore 
2700  Merced 
San  Leandro,  CA  94577 

1 Rockwell  International 
Autonetics  Missile  Sys  Div 
ATTN:  Dr.  M.  Chawla 
4300  E.  5th  Avenue 
Colunbus,  OH  43216 

3 , Sandia  Laboratories 

ATTN:  Dr.  W.  Herrmann 
„ Dr.  L.  Bertholf 
Dr.  J.  W.  Nunziato 
Albuquerque,  NM  87115 

1 Science  Applications,  Inc. 

101  Continental  Boulevard 
Suite  310 

El  Segundo,  CA  90245 

1 Science  Applications,  Inc. 

201  W.  Dyer  Road  (Unit  B) 
ATTN:  jG.  Burghart  ^ 

Santa  Ana,  CA  92707 
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5 Systems,  Science  § Software 
ATTN:  Dr.  R.  Sedgwick 
Ms . L . Hageman 
R.  G.  Herrmann 
D.  Grine 
M.  Rice 
P.  0.  Box  1620 
La  Jolla,  CA  92038 

I US  Steel  Corporation 
Research  Center 
125  Jamison  Lane 
Monroeville,  PA  15146 

1 Drexel  University 

Dept  of  Mechanical  Engineering 
ATTN : Dr . P . C . Chou 
32d  and  Chestnut  Streets 
Philadelphia,  PA  19104 

1 Forrestal  Research  Center 
Aeronautical  Engineering  Lab 
Princeton  University 
ATTN:  Dr.  A.  Eringen 
Princeton,  NJ  08540 

3 Southwest  Research  Institute 
Dept  of  Mechanical  Sciences 
ATTN:  Dr.  U.  Lindholm 
Dr.  W.  baker 
Dr.  P.  H.  Francis 
8500  Culebra  Road 
San  Antonio,  TX  78228 

3 Stanford  Research  Institute 
ATTN:  Dr.  L.  Seaman 
Dr.  D.  Curran 
Dr.  D.  Shockey 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


2 University  of  Arizona 
Civil  Engineering  Dept 
ATTN:  Dr.  A.  DaDeppo 
Dr.  R.  Richard 
Tucson,  AZ  85721 

4 University  of  California 
Los  Alamos  Scientific  Lab 
ATTN:  Dr.  R.  Karpp 
Dr.  J.  Dienes 
Dr.  L.  Germain 
Dr.  B.  Germain 
P.  0.  Box  808 
Livermore,  CA  94550 

1 University  of  Dayton 

Univ  of  Dayton  Research  Inst. 
ATTN:  Mr.  H.  F.  Swift 
Dayton,  OH  45405 

2 University  of  Delaware 

Dept  of  Mechanical  Engineering 
ATTN:  Prof.  J.  Vinson 

Dean  I . Greenfield 
Newark,  DE  19711 

1 University  of  Denver 
Denver  Research  Institute 
ATTN:  Mr.  R.  F.  Recht 

2390  S.  University  Boulevard 
Denver,  CO  80210 

2 University  of  Florida 
Department  of  Engineering 

Science  and  Mechanics 
ATTN:  Dr.  C.  A.  Sciammarella 
Dr.  L.  Malvern 
Gainesville,  FL  32601 

Aberdeen  Proving  Ground 

Cdr,  USATECOM 

ATTN:  Mr.  W.  Pless 

Mr.  S.  Keithley 
DRSTE-SG-H 
Dir,  USAMSAA 
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